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INTRODUCTION 


In  previous  reports  [1,  2,  3]  the  model  of  a  coupled  thermo-mechanical  be¬ 
haviour  of  ablative  composite  materials  under  intensive  heat  actions  was  devel¬ 
oped.  Different  aspects  of  this  model  in  detail  were  also  presented  in  our  works 
[5  -  10].  The  aim  of  the  previous  Report  was  to  verify  the  model  for  the  example 
of  Generic  Carbon-Phenolic  Textile  Fabric  2-D  Composite.  Experimental  results 
on  thermal  and  mechanical  properties  of  Carbon-Phenolic  Composite  were  given 
by  Wright  Lab.  [11].  As  the  result  of  the  comparison  conducted  of  the  theoretical 
and  experimental  results,  it  was  established  that  the  model  describes  a  behaviour 
of  ablative  composites  with  an  accuracy  being  sufficiently  acceptable  for  such  com¬ 
putations.  The  best  agreement  of  the  results  (error  w  7  - 10  %)  has  been  obtained 
for  heat-physical  properties: 

•  heat- conductivity, 

•  heat-capacity, 

•  gas-permeability,  and  also  for  mechanical  properties  (error  10  -  20  %)  at 
room  (294  K),  high  (600  -  1400  K)  and  superhigh  (1400  -  3000  K)  tempera¬ 
tures.  The  worse  accuracy  was  obtained  for  elevated  temperatures  (294  -  600 
K).  The  main  cause  of  this  fact  is  the  absence  of  reliable  experimental  data 
on  mechanical  properties  of  phenolic  matrix  under  high  temperatures,  as  in 
calculations  we  used  data  from  scientific  literature  and  also  we  did  not  take 
account  of  gradientness  (nonuniformity)  of  heating  under  these  temperatures 
although  the  heating  rate  was  sufficiently  high. 

The  aim  of  the  present  Report  is  to  conduct  the  next  stage  of  verification 
of  the  model:  prediction  of  a  behaviour  of  Generic  Carbon/Phenolic  Composite 
as  a  part  of  a  cylindrical  element  of  a  structure  under  non-uniform  heating.  A 
geometric  scheme  of  the  structure  element,  which  was  calculated,  was  proposed  by 
Wright  Lab.  Its  peculiarity  consists  in  that  principal  axes  of  orthotropy  of  Textile 
Composite  Material  contained  in  a  cylindrical  shell  are  not  coincident  with  axes  of 
the  cylindrical  coordinate  system  and  are  rotated  with  angle  7r/6  from  the  cylinder 
axis. 

It  is  evident,  that  thermo-mechanical  characteristics  of  such  turned  composite 
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(further  it  is  called  Carbon-Phenolic  Bias-Reinforced  Fabric  Composite  (BRFC)) 
in  the  cylindrical  coordinate  system  differ  from  the  corresponding  characteristics 
in  the  principal  axes  of  orthotropy.  Therefore,  the  purpose  of  Chapters  VII  and 
VIII  of  the  present  Report  (Chapters  I  -  VI  were  in  the  previous  report  [1])  is 
to  determine  heat-physical  and  mechanical  properties  of  BRFC  in  the  cylindrical 
coordinate  system. 

Chapter  IX  is  devoted  to  a  description  of  methods  for  calculation  of  an  abla¬ 
tion  rate  of  BRFC  under  intensive  non-imiform  heating.  A  main, attention  in  these 
calculations  is  paid  to  establishing  a  relationship  between  thermo-mechanical  and 
thermo-chemical  partials  within  a  whole  mechanism  of  ablation  of  BRFC. 

Calculation  of  the  ablation  rate  of  Carbon-Phenolic  Composite  according  to 
thermo-mechanical  mechanism,  mainly  for  sublimation  and  combustion,  was  in¬ 
vestigated  sufficiently  in  detail  up  to  now  (see,  for  example  [13]);  although  many 
important  questions  remain  still  uncleared  up.  For  example,  the  ablation  rate 
of  composite  materials  under  combined  diffusive-sublimational  regime  is  usually 
determined  only  in  a  phenomenological  way;  on  the  base  of  working  of  previous 
investigations. 

Efforts  in  prediction  of  thermo-mechanical  ablation  (mainly  mechanical  ero¬ 
sion  of  material  within  a  flow)  were  started  long  ago  [18].  However,  imperfection  of 
the  models  suggeted  before  (which  considered,  in  the  main,  only  thermo-stresses 
and  thermocracking  at  the  composite  surface)  did  not  allow  to  obtain  effective 
mathematical  methods  having  a  sufficient  accuracy  for  engineering  calculations. 

In  our  works  [3,  5]  the  model  has  been  suggested  for  calculation  of  thermo¬ 
mechanical  ablation  on  the  base  of  a  combined  analysis  of  thermo-mechanical 
stresses  and  physico-chemical  processes  occuring  in  composites  at  interaction  with 
high-temperature  chemically-reacting  gas  flows. 

In  §  9.1  of  Chapter  IX  there  are  results  of  calculations  of  the  ablation  rate  of 
BRFC  with  taking  account  of  both  the  mechanisms:  thermo-chemical  and  thermo¬ 
mechanical.  Comparison  of  a  relation  between  these  two  partials  has  been  con¬ 
ducted,  and  analysis  of  the  effect  of  Bias  angle  on  ablation  rate  of  BRFC  was  also 
performed. 

Chapter  X  is  devoted  to  calculation  of  heat  stresses  and  a  temperature  field 
within  a  cylindrical  shell  made  of  BRFC. 


Chapter  VII 

HEAT-PHYSICAL  PROPERTIES  OF 
BIAS-REINFORCED  FABRIC  COMPOSITE 

§7.1.  MODEL  OF  BIAS-REINFORCED  FABRIC 
COMPOSITE  (BRFC) 

Before,  in  previous  work  [1]  we  gave  thermo-mechanical  relations  of  the  model 
of  ablative  fabric  composite  in  principal  axes  of  orthotropy.  Fabric  composite,  in 
accordance  with  the  experimental  data  presented  by  Wright  Lab.,  is  assumed  to 
be  transversally  isotropic,  as  properties  of  the  composite  are  the  same  along  both 
the  warp  and  weft. 

In  cylindrical  elements  of  structures,  according  to  data  of  Wright  Lab.,  the 
fabric  composite  is  located  so  that  the  principal  axes  of  orthotropy  of  the  material 
Ox\  i  =  1, 2, 3  are  not  coincident  with  axes  of  the  cylindrical  coordinate  system  of 
the  structure:  namely,  the  axis  Ox'^  laying  in  the  direction  of  the  warp  is  inclined 
by  angle  ?/>  from  the  axis  Ox^  of  the  cylinder  (Fig.  1). 

Since  calculation  of  thermo-mechanical  processes  within  a  cylindrical  struc¬ 
ture  under  the  action  of  high-speed  hot  gas  flows  is  conducted  in  the  Oxi  axes  of 
the  cylindrical  coordinate  system,  it  is  necessary  to  obtain  design  characteristics 
of  the  fabric  composite  in  this  coordinate  system  Oxi. 

Further,  such  fabric  composite  with  a  inclined  direction  of  fabric  laying  is 
called  Bias-Reinforced  Fabric  Composite  (BRFC). 


§7.2.  THEORETICAL  RELATIONS  FOR 

HEAT-PHYSICAL  CONSTANTS  OF  BRFC 

Heat-  Conductivity. 

Let  be  components  of  a  heat-conductivity  tensor  of  BRFC  in  the  cylin¬ 
drical  coordinate  system  Oxi  i  =  1,2,3.  These  components  kijjp  are  connected  to 


A  scheme  of  orientation  of  principal  axes  of  orthotropy 
of  fabric  composite  in  a  cylindrical  structure. 


-  s- 


components  of  the  heat- conductivity  tensor  kij  in  the  principal  axes  of  orthotropy 
Ox'^  with  the  help  of  matrix  Qu'  of  rotation  about  the  0x2  axis  by  angle  ip: 


—  k%i  j'  Qii'Q  j'  j  5 

where  tensors  ki>ji,  kij^  and  Qn'  have  the  following  structure: 

(ki  0  0 

ki' j!  =  I  0  ^2  0 

\  0  0  fcs 

(kwxjj  0 

0  ^22^  ^ 

kl^tl;  0  ^33^ 


(1) 


(2) 


(3) 


(cos  Ip  0  —sin  Ip  \ 

0  10  (4) 

sin  Ip  0  cos  Ip  j 

Principal  values  of  the  heat-conductivity  tensor  k\ ,  k2  and  k^  are  expressed 
in  terms  of  characteristics  of  matrix  and  fibres  of  the  composite  by  formulae  (35) 
-  (38)  of  the  previous  Report  [1],  where  ^2  =  ^3  non-zero  values  of  the  heat- 
conductivity  tensor  of  BRFC  kn^,  A:22V>i  ^33V>  ki^^  axe  connected  to  ka  by 
the  relations: 

kill/)  =  ki  cos^  Ip  +  k2  sin^  ip] 

kzzrj,  =  ki  sin^  V’  +  ^2  cos^  ip]  (1') 


fci3V>  =  ^  sin2V’; 

k22tp  —  k2]  k 


h 


Heat  Capacity,  Density  and  Gas  Permeability. 

Heat  capacity  c  and  density  p  of  composite  BRFC  axe  scalar  characteristics 
and  independent  of  rotation  of  coordinate  axes,  and,  hence,  are  expressed  by 
formulae  (40)  and  (42)  from  [1]. 

Gas-permeability  tensor  Kij  of  fabric  composite  in  Ox\  axes  depend,  generally 
speaking,  on  orientation  of  coordinate  axes.  However,  in  our  previous  work  [1]  it 
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was  established  that,  in  fact,  gas  permeability  of  fabric  carbon-phenolic  composite 
is  the  same  in  the  Across  Ply  direction  K\  and  directions  of  the  warp  and  weft 
K2  =  A's,  i.e.  tensor  Kij  is  isotropic:  Kij  =  KiSij  (Ki  =  K2  =  K3). 

This  conclusion  is  also  verified  by  experimental  data  of  Wright  Lab.  [11]. 
Therefore,  in  rotation  of  the  coordinate  system  Oxi  tensor  Kij  does  not  change 
and,  hence,  gas-permeability  tensor  of  composite  BRFC  Kij^  is  coincident  with 
K^j: 

Kijjp  —  Kij  =  KiSij.  (5) 

Expression  for  A'l  was  given  in  [1]  (see  formula  (43)). 


Thermal  Expansion. 

Components  of  heat-expansion  tensor  of  composite  BRFC  in  Oxi  axes 
are  connected  to  components  Cij  of  heat-expansion  tensor  of  the  fabric  composite 
in  principal  axes  of  orthotropy  Ox'-  by  formulae  being  similar  to  (1)  and  (!’): 


Q  O 

S-ijip  —  S-i' j' Qii'Q j' j ■) 


(6) 


or: 


S-Wip  =  £1  cos^  V*  +  £3  sin^  V’) 
£33  v-  =  h  sin^  V’  +  £3  cos^  tp] 

/  o  o' 

£i3v>  =  I  — ;; —  I  sm2V>; 


(7) 


o  o  ,00, 

£22^”  —  ^3;  (£2  =  £3); 

where  tensors  Sij^  and  have  the  following  form: 

(£11^  0  £13^ 

o 

0  £22iI>  0 

o  o 

£i3V'  0  £22^ 


(8) 
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Expressions  for  e\  and  £2  were  given  in  [1]  (see  formulae  (28),  (29). 


§7.3.  CALCULATED  RESULTS 

In  order  to  analyze  the  influence  of  inclination  angle  V’  on  a  temperature 
dependence  of  heat  conductivities  of  the  composite  BRFC  fcii^(0)  and 
we  conducted  numerical  calculations  by  formulae  (!’)  for  different  values  of  angle 
=  0°,  15°,  30°,  15°,  60°.  Figure  2  shows  dependences  of  heat  conductivity 
of  BRFC  in  the  transverse  direction  Oxi  and  Figure  3  -  dependences  of  heat 
conductivity  kzzii>  in  the  longitudinal  direction  Oxz  on  temperature  9  of  heating 
for  different  values  of  the  angle  xj).  Here  and  further  the  heating  rate  is  assumed 
to  be  a  constant:  ^0  =  5  K/s.  Tables  1  and  2  present  corresponding  numerical 
values  of  and  kzzrp. 

With  increasing  angle  xp,  the  heat  conductivity  kn^  grows  but  kzztp  decreases. 
For  the  main  design  value  xp  =  30°  a  ratio  of  heat  conductivities  kn^p  and  ki  at 
xp  =  0°  is  equal  to  ku^^/ki  «  2  within  the  temperature  9  interval  from  300  K  to 
3000  K. 

Changing  angle  xp  has  a  less  effect  on  the  longitudinal  coefficient  of  heat 
conductivity  A;33^;  for  example  for  xp  =  30°,  the  ratio  kzz^/kz  is  equal  to  0.85  -  0.9 
within  all  the  temperature  interval. 

Analogous  investigations  on  the  effect  of  angle  xp  were  conducted  also  for  heat 
expansion.  Figures  4  and  5  show  results  of  numerical  computations  for  longitudinal 
£33^  and  transverse  £11^  heat  deformations  of  BRFC  depending  on  temperature 
9  of  heating  for  different  values  of  angle  xp  =  15°,  30°,  45°. 
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Pig.  2. 

Varying  Across  Ply  Thermal  Conductivity 
of  Carbon/Phenolic  BRFC  with  different  bias  angles  V’ 
versus  temperature  9  of  heating  at  5  K/s  rate. 


Across  Ply  Thermal  Conductivity  knf 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angles  'f' 


Temperature, 

9 

K 

Conductivity  W/(m  K) 

angle  ^ = 

15“ 

f=  30° 

45° 

340 

.775 

1.104 

1.553 

400 

.  779 

1.109 

1.560 

450 

.786 

1.117 

1.568 

500 

.805 

1.136 

1.588 

560 

.841 

1.171 

1.621 

615 

.867 

1.197 

1.647 

670 

.891 

1.223 

1.676 

800 

1.031 

1.385 

1.867 

1000 

1.675 

2.285 

3.119 

1200 

2.079 

3.029 

4.326 

1400 

2.401 

3.503 

5.008 

1600 

2.  794 

4.077 

5.829 

1800 

3.260 

4.  752 

6.  791 

2000 

3.769 

5.504 

7.874 

2200 

4.318 

6.328 

9.073 

2400 

4.916 

7.232 

10.395 

2600 

5.560 

8.214 

11.839 

2800 

6.246 

9.271 

13.402  ' 

3000 

6.973 

10.399 

15.080 
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Fig.  3. 

Varying  Longitudinal  Thermal  Conductivity  ^33^ 
of  Carbon/Phenolic  BRFC  with  different  bias  angles 
versus  temperature  0  of  heating  at  5  K/s  rate. 
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Tabl . 2 

Long’titudional  Thazmal  Conductivity 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angle  f 


Temperature,  K  Conductivity  W/(m  K) 

Q  _ _ ^ - 


Angle  (f' =15' 


340 

2.331 

395 

2.342 

450 

2.351 

505 

2.372 

560 

2.400 

615 

2.427 

670 

2.460 

800 

2.703 

1000 

4.563 

1200 

6.573 

1400 

7.615 

1600 

8.864 

1800 

10.323 

2000 

11.979 

2200 

13.828 

2400 

15. 876 

2600 

18.119 

2800 

20.557 

3000 

23.188 

f=  so""  45" 


2.002 

1.553 

2.012 

1.561 

2.020 

1.568 

2.040 

1.588 

2.071 

1.621 

2.097 

1.647 

2.129 

1.676 

2.349 

1.867 

3.953 

3.119 

5.  624 

4.326 

6.513 

5.008 

7.581 

5.829 

8.831 

6.  792 

10.244 

7.  874 

11.819 

9.073 

13.560 

10.396 

15.465 

11.840 

17.533 

13.402 

19.761 

15.081 
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Fig.  4. 

Varying  Across  Ply  Thermal  Expansion 
of  Carbon/Phenolic  BRFC  with  different  bias  angles  V’ 
versus  temperature  B  of  heating  at  5  K/s  rate. 


Tabl.  3 

o 

Across  Ply  Thermal  Expansion 
of  Carbon/ Phenolic  BRCM  with 
different  Bias  angle  if/ 


Temperature,  K  Expansion,  % 


P 


angle  15^ 

JO** 

45 

295. 

.045 

.037 

.025 

320. 

.270 

.219 

.151 

360. 

.628 

.511 

.350 

400. 

.983 

.  798 

.546 

440. 

1.312 

1.064 

.  726 

480. 

1.533 

1.241 

.844 

490. 

1.551 

1.256 

.853 

500. 

1.547 

1.252 

.849 

510. 

1.521 

1.230 

.833 

520. 

1.472 

1.190 

.805 

540. 

1.325 

1.069 

.  721 

560. 

1.184 

.954 

.640 

580. 

1.140 

.918 

.614 

600. 

1.213 

.976 

.653 

620. 

1.347 

1.085 

.  726 

640. 

.000 

.000 

.000 

660. 

1.627 

1.311 

.879 

680. 

1.749 

1.409 

.944 

700. 

1.850 

1.490 

.999 

720. 

1.926 

1.551 

1.040 

730. 

1.952 

1.572 

1.053 

740. 

1.969 

1.586 

1.062 

750. 

1.976 

1.592 

1.066 

760. 

1.973 

1.588 

1.063 

780. 

1.928 

1.552 

1.037 

800. 

1.829 

1.471 

.982 

820. 

1.450 

1.163 

.771 

840. 

.840 

.668 

.434 

880. 

.000 

.000 

.000 

920. 

.088 

.058 

.017 

960. 

-.638 

-.532 

-.388 

1000. 

-1.170 

-.966 

-.687 

1100. 

-1.479 

-1.219 

-.864 

1200. 

-1.346 

-1.112 

-.791 

1400. 

-1.842 

-1.510 

-1.056 

1600. 

-4.427 

-3.593 

-2.454 

1800. 

-8.481 

-6.862 

-4.650 

2000. 

-10.280 

-8.309 

-5.617 

2500. 

-10.212 

-8.253 

-5.576 

3000. 

-9.958 

-8.048 

-5.440 

temperature  0  ,  K 


Fig.  5. 

Varying  Longitudinal  Thermal  Expansion 
of  Carbon/Phenolic  BRFC  with  different  bias  angles  0 
versus  temperature  9  of  heating  at  5  K/s  rate. 


Tabl.  4 


Longitudional  Thermal  Expansion 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angles 


Temperature  f 

Q 

K 

Expansion^  % 

angle  ^=15° 

f  =  30“ 

O 

II 

295. 

.005 

.  014 

.025 

320. 

.031 

.  081 

.150 

360. 

.  069 

.187 

.348 

400. 

.103 

.289 

.543 

440. 

.132 

.381 

.722 

480. 

.145 

.438 

.839 

490. 

.144 

.441 

.847 

500. 

.141 

.438 

.844 

510. 

.135 

.427 

.827 

520. 

.126 

.410 

.  798 

540. 

.104 

.362 

.  714 

560. 

.  084 

.316 

.633 

580. 

.075 

.299 

.606 

600. 

.079 

.318 

.643 

620. 

.090 

.354 

.715 

640. 

.  000 

.  000 

.000 

660. 

.113 

.429 

.860 

680. 

.122 

.459 

.920 

700. 

.129 

.483 

.967 

720. 

.133 

.498 

.996 

730. 

.135 

.502 

1.003 

740. 

.135 

.502 

1.003 

750. 

.134 

.499 

.997 

760. 

.132 

.492 

.983 

780. 

.125 

.465 

.930 

800. 

.113 

.420 

.839 

820. 

.071 

.265 

.529 

840. 

.006 

.  026 

.052 

880. 

.000 

.  000 

.000 

920. 

-.075 

-.269 

-.535 

960. 

-.156 

-.563 

-1.119 

1000. 

-.217 

-.791 

-1.575 

1100. 

-.252 

-.978 

-1.971 

1200. 

-.233 

-.966 

-1.968 

1400. 

-.234 

-1.039 

-2.137 

1600. 

-.392 

-1.547 

-3.123 

1800. 

-.697 

-2.457 

-4.861 

2000. 

-.809 

-2. 830 

-5.592 

2500. 

-.706 

-2.711 

-5.450 

3000. 

-.594 

-2.568 

-5.264 
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Chapter  VIII 

MECHANICAL  PROPERTIES  OF 
BIAS-REINFORCED  FABRIC  COMPOSITES 


§8.1.  ELASTIC  MODULI 


Let  us  denote  now  Jijki  components  of  a  tensor  of  elastic  pliabilities  of  a 
fabric  composite  in  principal  axes  of  orthotropy  Ox[  and  Jijkiijj  in  turned  axes  of 
orthotropy  Oxi.  These  components  are  connected  by  the  relation  of  rotation: 

Jijklik  =  Ji'j'k'vQii'Qjj'Qu'Qkk'-  (9) 


From  the  relations  (9)  we  can  obtain  the  relations  connecting  elastic  moduli  of  fab¬ 
ric  composite  Ei,  Uij,  Gij  to  elastic  moduli  of  Bias-Reinforced  Fabric  Composites 

Eifp^  Gij'fp* 


Eii/j 


/  cos^  ^  sin^V’  sin^  2^^  f  ^ 


^  _  sin'^  ■0  cos^  0  sin^  20  i/i3  1 

E2jP  =  E2] 


•  o  /  /^H-2i^i3  ,  1  ^  ,  cos20 
+  _  1  +  _ 


Expression  for  elastic  moduli  Ei,  Uij,  Gij  depending  on  heating  temperature  9 
were  given  in  [1]  (relations  (24)). 

Figure  6  shows  calculated  data  of  a  dependence  of  elastic  modulus  E3^{$)  on 
heating  temperature  9  for  different  values  of  angle  0  =  15,  30  and  45°.  Table 
5  gives  corresponding  numerical  values  of  modulus  As  seen  from  the  Table 
5,  with  increasing  angle  0,  longitudinal  elasticity  modulus  Ez^  of  RRFC  sharply 
decreases:  for  angle  0  =  30°  modulus  Ez^,  is  equal  only  to  7  %  of  the  value  Ez  at 
9  =  295  K  and  0.6  %  at  2200  K.  It  can  be  explained  by  that  for  angles  0  >  5° 
elastic  properties  of  a  composite  are  defined,  in  the  main,  by  properties  of  its 
matrix. 


0  —  ■!  I  I  I  -J - 1 

300  600  900  1200  1500  1800  2100  2400  2700  3000 

temperature  6  ,  K 


Fig.  6. 

Varying  Longitudinal  Elastic  Modulus 
of  Carbon/Phenolic  BRFC  with  different  bias  angles  ^ 
versus  temperature  9  of  heating  at  5  K/s  rate. 


Tabl.  5 


Longitudional  Elastic  Modulus  £3,^ 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angles  V 


Temperature,  K  Modulus,  GPa 

9  _ 


angle  Y=0° 

15“ 

30“ 

Yk  45 

295.0 

60.228 

11.897 

4.437 

3.157 

670.0 

43.444 

3.464 

1.181 

.825 

920.0 

37.  615 

1.987 

.  664 

.461 

1365.0 

40. 691 

1.692 

.560 

.387 

2200.0 

35. 653 

.  708 

.229 

.156 

300  600  900  1200  1500  1800  2100  2400  2700  3000 


temperature  6  ,  K 


Fig.  7. 


Varying  Across  Ply  Elastic  Modulus  Ei^ 
of  Carbon/Phenolic  BRFC  with  different  bias  angles  'ip 
versus  temperature  6  of  heating  at  5  K/s  rate. 


Tabl.  6 


Across  Ply  Elastic  Modulus 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angles  yr 


angles:  0  ^^15"  Y=45‘‘  ^=75^ 


294.0 

6.  79 

6.14 

5.34 

5.59 

8.20 

20.16 

673.0 

1.75 

1. 59 

1.39 

1.47 

2.23 

6.42 

922.  0 

.96 

.  87 

.  77 

.82 

1.26 

3.77 

1366.0 

.  78 

.  72 

.64 

.69 

1.06 

3.24 

2200.  0 

.27 

.26 

.24 

.27 

.43 

1.38 

300  600  900  1200  1500  1800  2100  2400  2700  3000 


temperature  0  ,  K 


Fig.  8. 


Varying  Shear  Modulus 

of  Carbon/Phenolic  BRFC  with  different  bias  angles  V’ 
versus  temperature  9  of  heating  at  5  K/s  rate. 


Tabl.  7 


Shear  Elastic  Modulus,  GPa 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angles 


Temperature ,  K  Modulus,  GPa 


9 


angle  0* 

0 

II 

45° 

o 

Vo 

II 

295.  0 

1.679 

2.  098 

4.189 

8.353 

4.189 

675.  0 

.435 

.546 

1.121 

2.368 

1.122 

920.0 

.243 

.305 

.626 

1.321 

.626 

1365.0 

.204 

.256 

.522 

1.086 

.522 

2200.0 

.  083 

.  103 

.200 

.381 

.200 

Figure  7  and  Table  6  show  results  calculated  by  formula  (10)  Across  Ply 
Modulus  depending  on  heating  temperature  d  for  angles  ?/>  =  0,  15,  30,  45, 

60  and  75°.  As  seen  from  Figure  7,  the  effect  of  a  matrix  on  value  of  is  also 
prevalent  for  angles  =  0 . . .  60°.  Value  for  angles  V’  from  this  interval  changes 
very  weakly,  and  only  for  >  60°  function  starts  to  grow  approaching 

Ez{e). 

There  is  another  situation  for  shear  modulus  G\zri>-  Value  Gi3^  reaches  its 
maximum  for  angle  tf)  =  45°  (8.35  GPa  at  295  K  and  0.38  GPa  at  ^  =  2200  K), 
that  is  higher  by  four  times  than  for  fabric  composite  at  V’  =  0°.  For  angle  ^|^  =  30° 
this  modulus  is  higher  approximately  by  two  times. 


§8.2.  STRENGTH 

8.2.1.  ACROSS  PLY  STRENGTH 


Strength  condition  for  fabric  2-D  composite  in  accordance  with  the  developed 
modulus  is  determined  by  relations  (14)  and  (15)  from  [1].  In  order  to  obtain  an 
expression  for  ultimate  stress  for  which  there  occurs  destruction  of  BFRC 

in  tension  along  the  Oxi  axis,  it  is  necessary  to  substitute  the  following  relation 
into  equation  (15): 

O'!!  =  O'llV)  cos^ 

cri3  =  -^sin2V’;  >  0;  (11) 

(Ti2  =  0. 


As  a  result,  we  obtain; 

<^11^  = 


{ cos^  ij)  5m^20\ 


(12) 


Having  substituted  expression  (30)  from  [1]  in  place  of  Si{d)  and  5i3(0)  in  equa¬ 
tion  (12),  as  a  result  we  obtain  a  dependence  of  strength  5"^^  in  tension  in  the  axis 
Ox\  direction  on  temperature  9  of  heating. 

When  we  derived  equation  (12),  it  was  assumed  that  a  level  of  pore  gas 
pressure  p  is  negligibly  small.  The  effect  of  p  will  be  investigated  in  detail  below. 


300  600  900  1200  1500  1800  2100  2400  2700  3000 


temperature  0  ,  K 


Fig.  9  . 

Varying  Strength  in  Tension  along  Ox^  axis 
for  Carbon/Phenolic  BRFC  with  different  bias  angles  V’ 
versus  temperature  6  of  heating  at  5  K/s  rate. 


Tabl . 


Varying  Strength  in  Tension  along  axis 

for  Carbon/ Phenolic  BRCM 

with  different  Bias  angles  vs  temperature 
of  heating 


Temperature,  K 

P 

Strength, 

GPa 

angle  Y=0'‘ 

Y=15'‘ 

Y=3d' 

T=4^ 

Y=6Cf 

-^=75* 

294.0 

11.58 

12.32 

14.92 

21.03 

36.21 

86.92 

533.  0 

4.  85 

5.15 

6.24 

8.  77 

15.04 

35.87 

922.  0 

1.08 

1.15 

1.41 

2.02 

3.61 

9.27 

1644. 0 

.  71 

.  75 

.92 

1.33 

2.42 

6.48 

3033.0 

.46 

.49 

.60 

.88 

1.64 

4.66 

-13- 


As  follows  from  formula  (12),  destruction  of  composite  BRFC  in  tension  along 
the  Ox\  cixis  occurs  due  to  combined  normal  separation  of  fabric  layers  from  each 
other  along  the  axis  Ox'-^  and  interlayer  shear  in  the  inclined  plane  {x\x'^). 

Figure  9  shows  calculated  results  of  a  dependence  of  strength  on  temper¬ 

ature  6  of  heating  for  angles  V’  =  0,  15,  30,  45,  60  and  75°.  with  increasing  angle 
V’,  strength  5^  grows.  However,  for  angles  ip  from  the  interval  from  0°  to  30°  this 
increase  is  equal  only  to  30  %.  Essential  growth  of  the  strength  (approximately  by 
three  times)  is  for  angles  ip  >  60°.  Temperature  dependences  S^^{0)  for  all  angles 
Ip  are  similar  to  the  temperature  dependence  of  matrix  strength  ai(^)  (formula 
(34)  from  [1]):  under  tempeartures  from  the  interval  600  -  900  K  there  occurs 
sharp  falling  down  the  strength  cth^,  after  that  its  value  is  almost  stabilized. 


8.2.2.  STRENGTH  IN  TENSION  ALONG  0x3  AXIS 

Expression  for  strength  <733^  of  BRFC  in  tension  along  the  OX3  axis  can  be 
obtained  after  substituting  the  relationships  between  aij  and  (733^: 


=  (733^  sin^^; 


<733  =  <733^  cos^  ip\ 

(13) 

<713  =  ^0-33^  sin 2t/) 

into  equation  (14)  from  [1]: 

<733^  = 

(14) 

9+ 

/sin^i/>  sm^2'ip\ 

~  sr  ^  4S.3 )  ■ 

(15) 

This  equation  defines  destruction  of  composite  BRFC  by  type  (C)  (delamination 
of  composite).  If  we  substitute  relations  (13)  into  the  destruction  condition  by 
type  (A)  (15)  from  [1]  (breakage  of  fabric  threads  laid  in  the  reinforcing  direction 
Ox'2  or  0x3),  we  obtain: 


^33^  — 


c(^). 


c(^)  _  *^3 

COS^  ^  * 


(16) 


-30' 
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Varying  Strength  in  Tension  along  Oxs  axis 
for  Carbon/Phenolic  BRFC  with  different  bias  angles  ^ 
versus  temperature  6  of  heating  at  5  K/s  rate. 
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Tabl.  .3^ 

Strength  in  Tension  along  Ox^  axis 
of  Carbon/Phenolic  BRCM  with 
different  Bias  angles 


Tempera t ure ,  K 

9 

Strength 

,  MPa 

angle  =  0 

r=  5" 

V^^IO^" 

T=  15^ 

45' 

295.  0 

399.7 

132.9 

66.  7 

44.7 

23.2 

16.4 

920.  0 

87.2 

12.4 

6.2 

4.2 

2.2 

1.5 

2475.0 

84.9 

5.  8 

2.9 

2.0 

1.0 

0.  7 

3033.0 

76.  8 

5.3 

2.  6 

1.8 

0.9 

0.6 

0  900  1200  tSOO  1800  2100  2400  2700 

temperature  0  ,  K 


Fig.  11. 

Varying  Strength  in  Shear  in  Oxix^  plane 
for  Carbon/Phenolic  BRFC  with  different  bias  angles  V’ 
versus  temperature  9  of  heating  at  5  K/s  rate. 
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Tabl.  10 

Varying  Shear  Strength  in  0x13  Plane  for 
Carbon/Phenolic  BRCM  with  different 
Bias  angles  7^  vs  temperature  of  heating 


Temperature,  K 

9 

strength,  MPa 

angle  Y =0^ 

V=15'' 

-r=3o° 

-r=45" 

y=  ffo*’ 

r=75“ 

294 

25.14 

18.11 

12.92 

11.58 

12.92 

18.10 

533 

10.33 

7.52 

5.40 

4.85 

5.40 

7.52 

922 

2.  83 

1.80 

1.22 

1.08 

1.22 

1.80 

1644 

2.  06 

1.21 

.79 

.71 

.79 

1.21 

3033 

1.60 

.82 

.52 

.46 

.52 

.82 

-3Jf- 


Value  of  calculated  from  expression  (15),  as  a  rule,  is  essentially  higher  than 
5^  determined  by  formula  (14);  >  5^,  therefore,  the  destruction  condition 

for  BRFC  by  type  (A)  is  not  usually  realized.  Further,  formula  (14)  for  5"^ 
is  assumed  a  main  expression  for  ultimate  strength  in  tension  in  the  Oxz  axis 
direction. 

Figure  10  shows  graphs  of  dependence  of  strength  S^^{6)  on  temperature  9  of 
heating  for  different  values  of  angle  V’  =  0,  5,  10,  15,  30,  45°,  calculated  by  formula 
(15).  Table  9  gives  corresponding  numerical  values  of  the  functions.  As  seen  from 
these  graphs,  strength  5^  sharply  falls  down  with  increasing  the  angle  For 
■0  =  5°  it  decreases  by  three  times  under  0  =  295  K  and  almost  by  15  times  under 
9  =  3033  K.  It  can  be  explained  by  that  for  the  angle  0  =  5°  there  occurs  a  change 
of  a  destruction  type  of  fabric  composite:  if  for  0  =  0°  the  composite  destructs 
due  to  breakage  of  fabric  fibres  and  is  described  by  functions  and  (type 
(C)).  Since  strength  at  breakage  across  the  fabric  layers  and  at  interlayer  shear 
5i3  is  defined,  in  the  main,  by  properties  of  the  matrix,  for  angles  0  >  5°  strength 
S'^^(9)  behaves  itself  similarly  to  strength  of  the  matrix,  i.e.  sharply  decreases  in 
heating  up  to  800  K,  and  then  is  almost  stabilized  in  heating  up  to  3033  K. 


8.2.3.  STRENGTH  IN  SHEAR 

In  order  to  determine  strength  of  BRFC  in  shear  in  the  plane  Oxixz  (in  turned 
axes  Oxi),  it  is  necessary  to  substitute  the  following  expressions  for  stresses  aij 
for  this  case  of  loading  into  the  strength  condition  (15): 


(Til  =  -CTutp  sin 20; 


0-33  =  (Ti3^sin20; 


(17) 


<^13  =  <^13rP  COS  20. 

As  a  result,  we  obtain  that  destruction  in  shear  occurs  when  (Tizip  reaches  the 
ultimate  value  Snip: 


/  sin^  20  cos^  20  \ 


=  Snip'i  Snip 


-1/2 


(18) 
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Temperature  dependence  of  functions  and  5’i3(0),  as  before,  is  determined 

by  formulae  (30)  from  [1]. 

Equation  (18)  means  that  destruction  of  BRFC  for  all  angles  ij)  occurs  due 
to  delamination  of  the  composite,  and  thus  its  strength  under  high  temperatures 
is  defined,  in  the  main,  by  strength  of  its  matrix.  This  conclusion  is  verified 
by  calculated  results  presented  in  Figure  11  and  Table  10.  Calculations  for  the 
function  5'i3^(^)  were  conducted  by  formula  (18)  for  angles  ^  =  0,  15,  30  and 
45°.  Minimal  values  of  strength  5i3^  are  realized  at  V’  =  45°:  at  =  295  K 
the  strength  decreases  by  two  times  as  compared  with  the  values  for  ^  =  0°, 
and  at  0  =  3033  K  -  almost  by  four  times.  Value  at  ip  =  30°  is  close  to 
the  corresponding  value  at  7p  =  45°.  For  angles  =  45...75°  function  is 
summetrical:  =  5i3^(V>  —  7r/4),  i.e.  with  growing  angle  tp,  value  of 

increases. 
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Chapter  IX 

RATE  OF  ABLATION  OF  CARBON/PHENOLIC  BRFC 

§9.1.  VOLUMETRIC AL  AND  SURFACE  ABLATION  OF  BRFC 

Under  the  action  of  high-speed  hot  gas  flow  on  a  surface  of  structures  of 
aerospace  vehicles  made  of  composite  materials  there  occurs  their  ablation,  i.e. 
a  loss  of  mass.  In  accordance  with  the  developed  classification  [4,  9],  ablation 
processes  of  composites  are  divided  into  volumetrical  ablation  (pyrolysis,  thermo¬ 
decomposition,  thermo-oxidative  decomposition)  and  Unear  (surface)  ablation  (re¬ 
cession  of  substance  from  a  surface  of  composite). 

At  volumetrical  ablation  a  composite  density  decreases:  p{t)  <  p(0),  but  the 
volume  remains  practically  a  constant  V  ~  const,  and  at  purely  linear  ablation 
the  density  is  not  practically  change  p  «  const,  but  there  occurs  a  decrease  of  the 
material  volume  V{t)  <  V’(O).  In  reality,  there  can  occur  both  volumetrical  and 
linear  ablation  processes  in  composites  simultaneously. 

Main  quantitative  parameters  describing  ablation  processes  axe:  ablation  rate 
(mass  or  linear,  or  both  of  them)  and  specific  heat  of  ablation. 

Volumetrical  ablation  of  2-D  composite  in  accordance  with  the  developed  model 
is  described  by  mass  rates  of  ablation:  J/  of  fibre  and  Ji  of  matrix,  expressions 
for  which  are  determined  by  formulae  of  the  Arrhenius  type  (45)  from  [1]  and  by 
heat  of  ablation  Ae°  and  Ae°  and  by  gasification  coefficients  of  fibre  F /  and  of 
matrix  Fj. 

Numerical  values  of  these  constants  and  functions  corresponding  to  Carbon/ 
Phenolic  Generic  Composite  are  given  in  [1]. 

Changing  reinforcing  angle  ^  of  composite  has  no  effect  on  characteristics  of 
volumetrical  ablation,  therefore  for  BRFC  volumetrical  ablation  is  described  by 
the  same  parameters. 

In  the  present  work  let  us  consider  in  detail  formulae  for  the  rate  of  Unear 
ablation  D  of  BRFC.  But  in  the  beginning  let  us  systematize  the  formulae  for 
calculation  of  linear  rate  of  ablation  of  2-D  Fabric  Composite.  • 
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9.2.  LINEAR  RECESSION  RATE  OF 


2-D  CARBON/PHENOLIC  COMPOSITE 

Let  us  denote  f{x,  t)  =  0  equation  of  a  composite  surface  of  ablation  S(t)  mov¬ 
ing  in  three-dimensional  space  with  Descartes  coordinates  Oxi  and  basis  vectors  e^, 
then  under  the  condition  on  sufficient  smoothness  of  the  surface,  function  f{x,  t) 
satisfies  the  following  equation  [7]: 

^  +  D-|V/|  +  i?i;-V/  =  0  fsE(t), 

where  V  is  the  nabla- operator  over  coordinates  x,  is  the  velocity  of  mechanical 
motion  of  material  points  of  the  compositeon  its  surface  S  caused  by  mechanical 
forces  or  heat  expansion,  and  D  is  the  normal  rate  of  own  motion  of  the  ablation 
surface  D  caused  by  irreversible  physico-chemical  transformations  and: 

Velocity  v-£  is  determined  from  the  solution  of  the  general  problem  on  mechanics 
of  motion  of  solid  deformable  medium  with  a  mobile  boundary,  in  particular,  at 
small  deformations  of  composite  materials,  numerical  value  of  vj^  can  be  neglected 
as  compared  with  D.  To  calculate  the  rate  D  there  exist  two  main  approaches: 
from  the  point  of  view  of  1)  molecular- kinetic  theory  of  gas-fluid  media  [13]  and 
2)  thermodynamical  theory  developed,  in  the  main,  in  our  works  [5,  19].  Further 
the  second  approach  will  be  applied. 

In  accordance  with  the  developed  model  of  surface  ablation  processes,  the  rate  of 
surface  ablation  D  is  assumed  to  consist  of  a  superposition  of  velocities  of  different 
independent  ablation  processes: 

D  =  Ds  -l-  Dc  +  Df  -f-  Dm,  (19) 

where 

Ds  is  the  linear  rate  of  evaporation  (sublimation)  of  the  composite, 

Dc  is  the  linear  rate  of  chemical  reaction  of  separate  phases  of  the  composite  with 
reacting  gas  flow  (a  main  kind  of  this  type  of  ablation  is  combustion  of  carbon 
phases  of  the  composite:  char  and  carbon  fibre  in  oxygen  contained  in  gas  flow). 
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Dm  is  the  rate  of  thermomechanical  erosion  of  the  composite  under  the  action  of 
high-speed  gas  flow, 

Df  is  the  linear  rate  of  melting  of  several  phases  of  the  composite,  for  exam¬ 
ple,  of  particles  of  oxide  ceramics  AI2O3,  Si02  etc.  (For  the  considered  case  of 
Carbon/Phenolic  Composite  it  is  assumed  that  Df=0.) 

For  the  thermodynamical  approach  considered  here,  expressions  for  rates  Ds, 
Dc,  Dm  are  based  on  applying  the  second  law  of  thermodynamics  and  the  as¬ 
sumption  on  the  absence  of  energy  dissipation  on  surfaces  of  physico-chemical 
transformations. 


§9.3.  COMBUSTION  RATE 

Combustion  process  of  composite  carbon/phenolic  material  in  a  hot  gas  flow 
containing  chemically  active  components  (for  examle,  air)  is  a  phenomenon  of 
heterogeneous  combustion,  i.e.  a  process  of  chemical  reaction  of  the  type: 

(solid)  -f  (gas)  ^  (gas).  (20) 

As  a  solid  component  of  carbon/phenolic  composite  (without  account  of  special 
agents)  we  have  carbon  C  being  a  main  part  of  chemical  composition  of  the  car¬ 
bonized  phenolic  matrix  and  also  carbon  fibres. 

Therefore,  further  let  us  consider  a  multiphase  system  consisting  of  N  elements 
(j  =  1,...  ,iV)  wherein  j  =  1  is  solid  carbon  C  and  the  rest  j  =  2,. . .  ,N  are 
mixtures  of  components  of  the  overrunning  gas  flow  and  combustion  products.  All 
combustion  reactions  are  assumed  to  occur  on  the  surface  S  of  the  composite. 

It  is  follows  from  the  assumptions  that: 

(a)  on  the  one  side  of  the  separation  surface  S  (let  it  be  S'*")  there  is  only  component 
j  =  1,  on  the  other  side  (S“)  there  are  all  the  rest  components  j  —  2,...  ,N; 

(b)  if  Yj  is  the  mass  concentration  of  a  component  in  the  mixture,  then 

on  the  surface  S"*":  Yi  =  l,Yj  =  0  j  =  2, . . .  ,N; 

N 

the  surface  E“:  Y\  =  0,  Yj  =  1. 

i=2 


on 


(21) 


Further  we  will  assume  that  in  the  heterogeneous  combustion  process  there 
occurs  n  chemical  reactions,  for  which  we  can  write  the  stoichiometric  relations: 

N 

fc  =  l,...,n;  (22) 

i=i 

where  Vkj  are  stoichiometric  coefficients,  Mj  is  the  molecular  mass  of  the  _;-th 
component. 

Changing  concentrations  Yj  of  components  in  a  chemical  reaction  are  deter¬ 
mined  by  the  balance  law  of  components’  mass: 

PwDc{0-Yi)  +  Ii=mi-,  (23) 


ptoDc{Yj -0)  + Ij  =mj,  j=2,...,iV; 


(24) 


where  Dc  is  the  linear  rate  of  combustion  of  the  composite,  pu,  is  the  density  of  the 
composite  on  the  combustion  surface:  p^  =  pi<pi  +P3<p3,  Ij  is  the  flow  of  diffusion 
of  the  ^-th  component  from  or  to  the  reaction  surface  S,  rhj  is  the  mass  rate  of 
formation  of  the  j-th  component. 

In  relations  (23),  (24)  we  take  account  of  assumptions  (21)  on  heterogenity  of 
the  conbustion  process.  From  the  same  conditions  (21)  it  follows  that  li  =  1  and 
/i  =  0,  as  there  does  not  occur  a  diffusion  of  a  solid  component,  then  from  relation 
(23)  we  obtain: 


Pw 


(26) 


This  formula  allows  us  to  determine  the  linear  rate  of  combustion  of  the  composite, 
if  the  mass  rate  mi  is  known. 

Mass  rates  of  formation  of  components  rhj  are  determined  as  follows: 


*"j*^*’  i  =  1’  •  •  •  ’  (26) 

i=i 


where  is  the  mass  rate  of  k-th  chemical  reaction  which  is  determined  by  the 
law  of  acting  masses: 


here  symbol  Hf-i  means  a  product  of  concentrations  Yj  of  components  taking 
part  in  the  fc-th  reaction,  Ea  is  the  activation  energy.  Here  the  law  of  substance 
conservation  is  satisfied  for  the  system  as  the  whole: 


rhj  =  0.  (28) 

i=i 

Diffusion  flows  of  the  j-th  component  from  or  to  the  reaction  surface  are  deter¬ 
mined  with  the  help  of  the  so-called  Reynolds  heat-mass-transfer  analogy  [13]: 


(Kj  -  y„),  j  =  (29) 

where  (a/cp)  is  the  heat  transfer  coefficient,  Yje  is  the  concentration  of  the  j-th 
component  at  the  external  surface  of  the  boundary  layer  of  the  overrimning  flow. 
With  using  (29),  equations  (24)  can  be  written  in  the  form: 


rhj  +  {a/cp)Yje 

PwJ^c  "H  is^l ^p) 


3=2,...  ,N. 


(30) 


The  equations  (30)  together  with  (26)  are  the  system  of  (2iV  —  1)  algebraic 
equations  to  determine  (2iV  —  1)  unknown  functions  rhi,  ihj,  Yj  j  =  2, ... 
on  solving  the  system  we  And  the  desired  linear  rate  of  combustion  Dc  by  formula 
(25). 

Let  us  consider  now  in  detail  a  chemical  composition  of  an  overrunning  flow. 
Usually  under  high  temperatures  (5000  K  and  higher),  the  overrunning  flow  disso¬ 
ciates:  besides  ’heavy’  components  O2,  N2  there  are  also  ’light’  ones  O,  N  therein. 
Let  us  assume  that  on  the  external  surface  of  the  boundary  layer  of  the  overrun¬ 
ning  flow  gas  consists  of  five  components:  its  composition  is  O,  O2  N,  N2,  NO 
(indexes  of  the  components  are  j  =  2, 3, 4, 5, 6).  On  the  composite  surface  S“  due 
to  combustion  reaction  the  gas  mixture  contains  also  components:  CO,  CO2,  CN 
and  C2N2  (i  =  7, 8, 9, 10).  At  the  combustion  surface  2  the  following  chamical 
reactions  run: 

^7-1-0  =  CO  reaction  =  1; 

C  +  O2  =  CO2  reaction  k  =  2; 

C  +  N  =  CN  reaction  fc  =  3;  (31) 


2C  +  N2  =  C2N2 


reaction  k  =  4. 


Chemical  composition  of  such  10-component  mixture  (component  j  =  1  is  car¬ 
bon  C  of  the  solid  phase)  and  chemical  reactions  (31)  allow  us  to  find  numerical 
values  of  molecular  masses  Mj ,  stoichiometric  coefficients  and  cilso  concentra¬ 
tion  values  Yje,  and: 

Yje=0  i  =  l,7...10.  (32) 

Number  n  for  such  mixture  is  equal  to  n  =  4.  Mass  rates  of  chemical  reactions  Uk 
(27)  are  written  as  follows: 


exp 


k  =  ,4, 


and  mass  rate  mi  of  vanishing  carbon  due  to  relations  (26)  and  (33)  takes  the 
form: 

m,  =  M,  ^  unUtY^.,  exp  ■  (34) 


Concentrations  Yk  contained  in  this  expression  can  be  determined  from  equation 
(30)  whereinto  we  should  substitute  relations  (25)  and  (33): 

^  Mi  EL,  exp  +  (a/c,)yj, 

Y^.  _  _ _ I _ V _ 

Efc=i  i'ikYk+iU^+^  exp  +  i<^/cp) 


i=2,...  ,iV. 


The  equation  (35)  is  the  system  of  nine  nonlinear  algebraic  equations  to  determine 
nine  concentrations  of  components  Yj{9y^)  in  the  form  of  functions  of  the  surface 
temperature 

Expression  for  the  linear  rate  of  combustion  Dc  as  a  function  of  is  determined 

by  formulae  (25)  and  (34)  after  finding  the  concentrations  Yj(d)  from  (35). 

In  accordance  with  the  method  described  dabove,  numerical  calculations  were 
conducted  for  combustion  rate  Dc  of  Carbon/Phenolic  Composite.  All  the  con¬ 
stants  Yje,  Mj,  Uk  and  Eak  contained  in  equations  (25),  (32)  -  (35)  were  taken 
from  data  [13]. 

Figure  12  shows  calculated  results  of  linear  rate  of  combustion  Dc  for  2-D 
Carbon/Phenolic  Composite  in  the  form  of  dependence  Dc{9w)  on  the  surface 


temperature  for  the  case  of  p  =  5  atm,  heat  transfer  coefficient  (a/cp)  was  chosen 
as  (a/cp)  =  2.2  kg/m^-  s.  As  seen  from  Fig.l2,  the  function  Dc(^u>)  has  a  typical 
character  for  combustion  processes  of  caxbon  materials  [20].  We  can  observe  two 
stages  of  combustion:  kinetic  stage  when  combustion  rate  Dc  sharply  increases 
(temperature  interval  is  from  1000  to  1700  K)  and  diffusive  stage  when  combustion 
rate  Dc  is  stabilized  at  the  certain  level  (temperature  interval  is  from  1700  to  3000 
K).  In  the  diffusive  regime,  the  combustion  rate  magnitude  is  determined,  in  the 
main,  by  rates  of  diffusion  of  oxidizers  O2  and  0  to  the  combustion  surface  S  and 
of  combustion  products  CO,  CO2  from  the  combustion  surface  S  of  the  composite. 

Under  higher  temperatures  d^u  of  heating  a  sublimational  regime  of  combustion 
of  the  composite  starts. 


§9.4.  SUBLIMATION  RATE 

Under  superhigh  temperatures  (higher  than  3000  K)  there  occurs  evaporation 
(sublimation)  of  solid  carbon  phases  of  the  composite:  carbon  fibres  and  car¬ 
bonized  matrix 

C(solid)  — >  C(gas).  (36) 

In  this  case  there  is  a  gaseous  component  j  =  1  at  the  combustion  surface  S 
and  the  equation  of  mass  balance  for  component  j  =  1  (of  the  type  (23))  can  be 
rewritten  as  follows: 

p^D,(Yi  -  1)  +  (^^^  (Yi  -  0)  =  0.  (37) 

here  Y\  is  the  concentration  of  gaseous  C  at  the  surface  S“;  and  at  the  external 
surface  of  the  boundary  layer  the  concentration  of  products  of  physico-chemical 
transformations,  as  before,  is  equal  to  zero;  Ds  is  the  linear  rate  of  sublimation. 
From  (37)  we  find: 

Y  _  _ Pw^s _ 

pwD s  ^p) 

The  rest  components  of  the  gas  mixture  j  =  2, ...  ,N  at  the  surface  S”  for 
the  given  case  can  be  considered  as  inert  with  respect  to  the  sublimation  process: 


mm/s 


temperature  0  ,  K 


Fig.  12. 

Varying  Rates  of  Thermo-chemical  (Sublimation  Dg 
and  Combustion  Dc)  Recession  of  2-D  Carbon/Phenolic  Composite 
versus  temperature  6  under  p  =  5.25  atm. 
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Relation  between  Different  Parts  of 
Thermal -Chemical  Ablation 
(Combustion  and  Sublimation) 
for  2-D  Carbon/Phenolic  Composite 
at  Pg-B  atm  in  air 


Tempera t ure  ,  K 

9 

Rate 

of  Recession,  mm/s 

Combustion 

2>c 

Sublimation 

2>s 

Total 

3>7C 

1500.0 

.000 

.047 

.047 

2000.0 

.000 

.127 

.127 

3000.0 

.010 

.136 

.146 

3200.0 

.052 

.136 

.188 

3400.0 

.224 

.136 

.361 

3600.0 

.828 

.136 

.964 

The  sublimation  process  can  be  considered  as  a  phase  transformation  of  the 
first  type.  For  phase  transformations  of  type  (36)  an  expression  for  the  phase 
transformation  rate  being,  in  fact,  a  linear  rate  of  sublimation  Ds  is  determined 
by  the  second  law  of  thermodynamics,  if  we  made  the  assumption  on  the  absence 
of  sources  of  energy  dissipation  on  the  surface  of  the  phase  transformation.  In  this 
case,  as  follows  from  [19],  Ds  is  determined  as: 

(39) 


PwDs  — 


/ 


V 


(p*5„  -cr„) 


V  P9 


\ 


1/2 


SnMj/ 


Sn  =  exp 


f  EasMA 

V  RaSr,  /’ 


(40) 


where  (7„  for  the  considered  case  is  the  partial  pressure  of  the  gas  phase  on  the 
surface  E: 


—  Pe^l ) 


(41) 


and  p*  is  the  constant  characterizing  the  pressure  of  ’initial  state’  of  the  gas  phase 
generated,  and  for  the  constant  the  following  relation: 


P*gSnMi 


<  1 


is  valid,  then  expression  (39)  takes  the  form: 


PwDs  —  ^PwP*g  ^Sn 

This  equation  is  similar  to  the  Lengmeur-Knudsen  equation. 

Two  equations  (38)  and  (42)  allow  us  to  evaluate  the  rate  D3  and  concentration 
Fi  in  the  form  of  functions  of  the  surface  temperature  Ojjj.  The  linear  rate  of 
sublimation  is  usually  relatively  small  so  that  the  relation  ^PgSnia/cpY /{p^pD  < 
1  is  satisfied,  then  equation  system  (42),  (38)  admits  an  approximate  solution  in 
the  form: 


temperature  K 


Fig.  13. 

Varying  Rates  of  Thermo-chemical  (Sublimation  Ds 
and  Combustion  Dc)  Recession  Dtc  of  2-D  Carbon/Phenolic  Composite 
versus  temperature  0^  '  '  and 

different  static  pressures:  1,  10,  100  and  500  atm. 
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Varying  Rates  of  Thermal -Chemical 
Ablation  of  2-D  Carbon /Phenolic 
Composite  on  Temperature  in  air, 
at  different  values  of  static  pressure 


Tempera ture , 

(AT» 


D- 


TC 


Rate  of  Recession,  mm/s 


P=1  atm  10  atm  100  atm 

e 

theory  experim.  theory  experim.  theory  experim. 


2000.0 

.13 

.15 

3000.0 

.14 

.15 

3200. 0 

.15 

.19 

3400.0 

.22 

.22 

3600.0 

.43 

.30 

3800.0 

1.07 

.  60 

4000.0 

2.  80 

3.  00 

4200.0 

4400.0 

4600.  0 
4800.  0 

5000.0 


.13 

.15 

.13 

.15 

.14 

.15 

.14 

.15 

.14 

.15 

.14 

.15 

.14 

.19 

.14 

.15 

.17 

.20 

.14 

.15 

.23 

.26 

.15 

.20 

.40 

.40 

.16 

.25 

.82 

.  72 

.20 

.30 

1.77 

2.10 

.30 

.40 

.50 

.60 

.88 

1.50 

1.58 

2.30 

As  follows  from  this  equation,  dependence  of  Dg  on  pressure  pe  of  tho  surroundings 
is  inversely  proportional.  This  expression  will  be  consider  further  as  a  sublima¬ 
tion  rate  for  Carbon/Phenolic  Composites.  Constants  p*  and  Eas  iii  (43)  are 
determined  in  experiments  and  are  equal  to: 


Eas 

RaM\ 


=  80  •  10^ 


K’ 


(44) 


Figme  12  shows  a  dependence  of  function  Ds{6w)  on  the  surface  temperature 
of  Carbon/Phenolic  Composite  determined  by  equation  (44).  Figure  12  also 
shows  a  dependence  of  summarized  thermochemical  rate  of  ablation  Dtc- 


Dtc  =  Ds  +  Dc 


(45) 


of  Carbon/Phenolic  Composite  on  temperature  0s  for  the  case  when  pe  =  Po  =  5 
atm.  As  seen  from  this  figure,  in  heating  up  to  3000  K  only  kinetic  and  diffusive 
regimes  of  combustion  of  the  composite  are  realized  in  air  and  Dc  ~  0,  and  imder 
higher  temperatures  (>  3000JC)  there  occurs  a  sublimation  regime  wherein  the 
ablation  rate  Dtc  sharply  increases  due  to  fast  growing  the  subhmation  rate  Dg . 

It  is  well  known  that  the  sublimation  rate  depends  on  total  pressure  pe  of  gas 
flow:  with  increasing  pressure  pe  the  rate  Ds  lowers.  Figure  13  gives  calculated 
results  of  the  rate  of  thermochemical  ablation  Dtc{^w)  for  three  values  of  pressure 
Pe  =  1,  10  and  100  atm  and  corresponding  experimental  data  from  work  [20]. 

Numerical  results  of  the  calculations  and  experiments  are  presented  in  Table 
14.  As  seen  from  the  table  and  Fig.  13,  formula  (44)  describes  the  dependence  of 
Ds  on  pressure  pe  not  only  qualitatively  but  allows  to  obtain  a  sufficient  accuracy 
of  quantitative  description. 


§9.5.  THERMO-MECHANICAL  ABLATION  RATE 
9.5.1.  GENERAL  RELATIONS 

Thermo-mechanical  ablation  means  the  process  of  mechanical  crumbling  off 
(erosion)  of  particles  from  a  surface  of  composite  by  a  hot  high-speed  flow. 
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The  cause  of  thermo-mechanical  ablation  can  be  main  acting  factors:  (1)  over¬ 
running  gas  flow  (with  velocity  V),  the  boundary  layer  of  which  moves  along  the 
ablation  surface  S  of  the  composite  with  velocity  V^\ 

Vnfi,  (46) 

where  =  V  •  n  is  the  normal  component  of  velocity  V,  and  (2)  internal  flow  of 
gaseous  products  of  pyrolysis  of  the  matrix  which  moves  through  pores  towards 
the  external  surface  2  of  the  composite  and  ’breaks  up’  a  thin  nonstrong  coke 
layer  of  the  composite  preventing  partially  from  gas  filtration.  Therefore,  the  rate 
of  thermo-mechanical  ablation  Dm  can  be  represented  as  a  sum  of  rates  and 
Di  of  these  two  processes: 

Dm  =  Dc  -t-  Di,  (47) 

where  Dg  will  be  called  a  rate  of  external  thermo-mechanical  erosion,  and  Di  -  & 
rate  of  internal  thermo-mechanical  erosion. 

From  the  definition  of  thermo-mechanical  erosion  it  follows  that  rates  Dg  and 
Di  should  depend  on  characteristics  of  resistance  to  the  force  action  of  external 
and  internal  gas  flows,  i.e.  on  strength  of  the  composite.  However,  the  composite 
strength  depends,  generally  speaking,  on  the  direction  of  mechanical  force  action, 
hence,  the  linear  rate  of  thermo-mechanical  ablation  also  depends  on  the  direction 
of  motion  of  a  gas  flow  with  respect  to  principal  orthotropy  axes  of  the  compos¬ 
ite.  It  is  a  main  distinction  of  rate  Dm  from  rate  of  thermo-chemical  ablation 
Dtc  which  is  independent  of  both  the  direction  of  flow  motion  and  orientation  of 
orthotropy  axes  of  the  composite. 

Thus,  linear  rate  of  thermo-mechanical  ablation  Dm  should  depend  on  some 
tensors  x  =  1, ...  ,n,  characterizing  a  type  of  anisotropy  of  a  composite 

material  (ositropic,  transversally  isotropic,  orthotropic  etc.),  and  on  the  orientation 
of  anisotropy  axes  with  respect  to  the  ablation  surface  2,  i.e.,  in  fact,  on  normal 
vector  n  to  the  surface: 


Dg=Dg(n,a^^\Vj:,e,,y, 


(48) 


Di  =  , 


X  =  1, . . .  ,  n. 


(49) 
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Besides  that,  rate  of  carring  off  Dt  depends,  evidently,  on  the  surface  temper¬ 
ature  9w  and  also  on  both  rate  and  direction  of  motion  of  hot  gas  flow  along 
the  ablation  surface  S  that  is  presented  in  formula  (48). 

Since  the  function  Dm  is  scalar,  arguments  in  (48),  (49)  should  also  be  scalar 
invariants  of  the  values  n,  Fs-  Let  us  consider  the  following  invariants: 


PS  =  Pge^S  '  Vt.  =  PgeVl\  >f  =  1, .  .  .  ,  n.  (50) 

Here  we  assumed  that  all  tensors  are  of  the  second  rank,  and  pg^  is  the 

gas  density  on  the  external  surface  of  the  boundary  layer  of  the  overrimning  flow 
(scalar). 

Number  of  invariants  n,  as  it  was  shown  in  [21],  is  not  more  than  6  (n  <  6). 


Hence,  rates  and  Di  can  be  presented  in  the  form  of  function  of  two  invariants 
and  p^^L 


De  —  De  \  a^^\p^,6wj  ; 

(48') 

Di  =  Di(a^^\erv);  X=l,...,n. 

(49') 

Now  let 

US  assume  that  in  (48’),  (49’)  the  dependence  on  arguments 

is  only 

linear: 

n 

De  =  Dex  {PE,9w) 

>f=l 

(51) 

n 

(52) 

>t=\ 


and  investigate  the  effect  of  the  type  of  anisotropy  on  the  rates  of  thermo-mechanicall 
ablation  Dg  and  Di. 


9.5.2.  ISOTROPIC  COMPOSITES 

For  isotropic  materials,  for  example,  nonreinforced  matrices, 'composites  re¬ 
inforced  by  dispersed  particles  or  reinforced  chaotically  by  short  fibres,  number 
n  =  1,  and  tensor  is  single  and  invariants  (50)  have  the  form: 

n  =  1,  =  ^,  a***  =  n  ■  6  ■  n  =  1; 


(53) 
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and  expressions  (51),  (52)  for  De  and  Di  take  the  form: 

De  —  DMe  —  elipT,  1  ^  ■w)‘i  (^4) 

Di  =  DMi  =  Dix{e^)-  (55) 

where  8  is  the  Kronecker  tensor  with  components  8ij  in  the  Oxi  axes.  Thus,  as  it 
should  be  expected,  for  isotropic  media,  rates  of  thermo-mechanical  ablation  Dm^. 
amd  jDmi  are  independent  of  direction  of  the  normal  n  and  are  determined  only 
by  the  surface  temperature  O^v  and  invariant  which,  in  fact,  is  a  velocity  head 
of  the  gas  flow  on  the  ablation  surface  S. 

Formulae  (54)  and  (55)  establish  only  on  what  arguments  the  rates  Dmc  and 
DMi  of  matrix  ablation  are  dependent.  Further  content  of  the  paragraph  is  devoted 
to  determination  of  the  functions  DMeiP'Li^w)  and  DMii&w)- 


RATE  OF  INTERNAL  EROSION  OF  MATRIX  Dni 


Expression  for  rate  Di  is  determined  on  the  basis  of  joint  solving  the  heat- 
cunductivity,  filtration  equations  and  kinetic  equation  describing  the  pyrolysis 
process  in  quasistatic  approximation  (^  -^  0)  and  has  the  form: 


dip  2 
dx 


Ji 

■p2exp{ 


pDi 


F^atn  \ 


0  <  X  <  -|-oo; 


—  _|.  d^$  _  ^ 

dx  PmCmDi  dx"^ 

dpg{\  -  P2)  K  d'^pg  dp2  ^ 

dx  Di  dx"^  ^  dx 


(56) 


Pg=Pe]  0  =  $^,  x  =  0;  0  =  Sq,  <P2  =  P2^  X  = -foo; 


dpgi'i-  -  p) 
dx 


where  Oq  and  p\  —  \  are  the  temperature  and  polymer  phase  concentration  of  the 
cold  material;  x  -  mobile  coordinate  connected  to  the  ablation  surface:  x  —  x  —  xd 


(Fig.  20a),  pg  -  density  of  the  gas  phase,  -  pre-exponential  multiplier,  K  -  gas 
permeabliity  coefficient  of  the  material  q;iatrix. 


Models  of  processes  of  thermo-mechanical  erosion  of  composites 
for  calculations  of  linear  ablation  rates: 

a)  Model  of  internal  erosion  of  a  matrix  for  calculation  of  rate  Dmi 

b)  Model  of  external  erosion  of  a  matrix  for  calculation  of  rate  Dmb 

1  -  matrix  (solid  phase),  2  -  pore,  3  -  dispresed  particles  of  matrix  in  gas  flow, 
4  -fibre,  5  -  pore  gas  pressure,  6  -  overrunning  hot  gas  flow. 


Models  of  processes  of  thermo-mechanical  erosion  of  composites 
for  calculations  of  linear  ablation  rates: 

c)  Model  of  external  erosion  of  unidirectio»al  composite  for  calculation  of  rate  D, 

d)  Model  of  external  erosion  of  unidirectional  composite  for  calculation  of  rate  D 

e)  Model  of  external  erosion  of  fabric  composite  for  calculation  of  rate  D^\ 

1  -  matrix  (solid  phase),  2  -  pore,  3  -  dispresed  particles  of  matrix  in  gas  flow, 


A  typical  graph  of  dependence  of  pore  gas  pressure  Pg  on  coordinate  x  for 
ablative  matrix  has  a  local  maximum  (Fig.  20a).  The  presence  of  this  maximum 
is  a  typical  peculiarity  of  heating  of  materials,  wherein  there  can  occur  a  thermo¬ 
decomposition,  that  was  established  by  direct  numerical  solving  the  equations  (56) 
of  internal  heat-mass-transfer  [1,6].  In  order  to  find  an  approximate  analytical 
solution  of  system  (56),  a  boundary  condition  for  the  filtration  equation  of  the 
problem  (56),  with  using  the  existence  of  local  extremum  of  function  p{x)  at  a 
certain  point  x*,  was  formulated  in  the  form:  dpfdx  =  0  at  a:  =  x*.  At  the  same 
point  X*  a  thermo- decomposition  process  only  starts,  therefore  we  can  assume  that 
<P2{X*)  1.  Solution  of  the  third  equation  of  system  (56)  is  sought  within  the 

domain  0  <  x  <  x*. 

Then  solution  of  the  system  (56)  has  the  form: 

^  =  ^(1 -^2)(Pm -/>s),  0<x<x*; 

V?2  =  exp >  0  <  X  <  -t-oo;  (57) 

e-9Q=  6w  exp(- , 

where 

/(X)  ^  exp  (58) 

Herein  it  is  taken  into  account  that  usually  pg  p  for  composite  materials  and 
1-(^2(X*)<1. 

To  determine  integral  I  the  following  approximation  is  used: 

- ) 


In  particular,  we  obtain  that 


T/n\  kjfiRa^w  f  Ram\ 

\  R»J' 


If  (^2  =  <Pw  denotes  a  value  of  solid  phase  concentration  of  the  matrix  at  the 
ablation  surface  x  =  0,  then  from  equation  (59)  at  y  =  0  we  obtain  the  expression 
for  combustion  rate  Dii 

D.  = 

Pm  In  (Pw 


(60) 
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In  accordance  with  the  considered  model,  ablative  material  at  its  external  sur¬ 
face  of  combustion  has  a  regular  structure  (Fig.  20a)  in  the  form  of  cells  with  wall 
thickness  hy,  and  cavity  dimension  Iq  connected  to  by  the  relation: 

(l  +  =  1  -  ifw  (61) 

Pyrolise  process  leads  to  descreasing  wall  thickness  hw  However,  the  value 
hyj  =  0  can  not  be  reached,  as  at  a  certain  value  of  hy,  there  occurs  mechanical 
destruction  of  the  walls  under  the  action  of  bending  stresses  a  in  the  wall  caused 
by  pore  pressure  Pg  of  pyrolytical  gases. 

An  idealized  scheme  of  a  stress  state  of  each  the  bulge  of  the  ablation  surface 
can  be  represented  as  cantilever  bending  of  a  beam  by  uniform  pressure  pg.  When 
the  bending  stress  exceeds  the  strength  limit  of  the  matrix  in  tension,  there  occurs 
breaking-off  a  surface  particle  and  its  carrying-off  to  the  gas  flow. 

The  considered  idealized  model  of  erosion  assumes  that  at  each  the  elementary 
destruction  a  shape  of  a  dispersed  particle  and  a  shape  of  the  residual  ablation 
surface  remain  self-similar  (Fig.  20a). 

The  value  of  for  which  there  occurs  destruction  of  the  next  layer  of  the 
matrix  surface  is  determined  by  the  formula  for  bending  stresses: 


hwy  _  6{pg  -pe) 
k  '  CTmT 


(62) 


where  cTmT  is  the  ultimate  stress  in  tension  (strength)  of  the  matrix  under  tem¬ 
perature  0yj. 

Expression  for  pressure  difference  (pg  —  Pe)  near  the  ablation  surface  can  be 
found  with  the  help  of  linearized  formula  (57): 


Pg  Pe  —  (^^) 

where  it  was  taken  into  account  that  at  the  combustion  surface  gas  density  in 
pores  pgy,  is  much  less  than  density  of  solid  phase  pgy,  <C  Pm- 

On  substituting  formula  (63)  into  (62)  and  then  (62)  into  (61),  the  equation  for 
Py^  can  be  derived: 


^DiPjnlo 

KffmT 


(64) 


As  a  rule,  (6Diplo)({Kao)  >  1,  therefore,  solution  of  equation  (64)  can  be  repre¬ 
sented  in  the  form: 

Since  ipw  <  1,  we  can  make  one  more  approximation,  namely 


KcTmT  \3/5 


(66) 


On  substituting  formula  (59)  and  (66)  into  (60),  the  final  expression  for  the 
erosion  rate  Di  can  be  found: 

2u;—l  /  r? 

/  JXa^xi 


n  —  T)  —  (  \  \ 


w  \  f  Ea^ 
exp 


Rad 


U3 


w , 


=  5/7. 
(67) 


RATE  OF  EXTERNAL  EROSION  OF  MATRIX  Dui 


Unlike  internal  erosion  when  material  particles  at  the  surface  destruct  under  the 
action  of  internal  forces  (intrapore  pressure),  ablation  surface  particles’  destruction 
(erosion)  at  the  presence  of  an  external  high-speed  flow  occurs  due  to  the  action 
of  local  pressure  ps  onto  bulging  parts  of  the  ablation  surface  (Fig.  20b). 

A  scheme  of  mechanical  destruction  of  an  matrix  surface  under  the  action  of 
an  external  flow  is  analogous  to  the  scheme  of  internal  erosion,  therefore  formulae 
(58)  -  (62)  also  take  place  but  instead  of  the  pressure  difference  Pg  -pe  we  should 
consider  the  pressure  head  ps : 

6pE 


^mT 


(68) 


On  substituting  formula  (68)  into  (62)  and  then  into  (61),  we  find  the  expression 
for  function  pw- 


-In  (pu,  =  X  = -ln(l  -  (1  +  ^), 

\crmT  /  ^  <^mT  ' 

Then  from  equations  (59),  (60)  and  (69)  we  obtain  the  expression  for  Bm^'- 


(69) 


But 


Pm 


\Eam  J 


1/2 


"  m  1 

iX(Psl<^mT)) 


Ea 


exp 


2Rey 


(70) 
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9.5.3.  TRANSVERSALLY  ISOTROPIC  COMPOSITES 

For  transversally  isotropic  media  which  have  only  one  axis,  for  example  OX3,  of 
transverse  isotropy  (they  are,  for  example,  unidirectionally  reinforced  composites), 
a  number  n  =  2,  there  are  two  tensors  and  invariants  (50)  have  the  form: 

n  =  2,  a\j^  =  SizSjz't  =  SnSji  +  ^t2<5j2; 

a(2)=nl,  a(i)=n?  +  nl,  (71) 

and  expressions  (51),  (52)  for  the  rate  De  take  the  form: 

De  =  ^||(pE,^«.)n|  +  Dx(pE,^u>)(nf  +nl).  (72) 

As  follows  from  formula  (72),  for  transversally  isotropic  composites  there  exist 
two  different  functions:  -  rate  of  external  erosion  of  surface  E  being  orthogonal 

to  the  0x3  axis  of  transverse  isotropy  and  Dx  -  rate  of  external  erosion  of  surface 
S  of  the  composite  being  orthogonal  to  axis  Oxi  or  0x2. 


RATE  OF  THERMO-MECHANICAL  EROSION  OF  FIBRES 

Although  fibres  without  matrix  are  not  independent  constructional  material, 
but  theoretically  we  can  consider  separately  a  process  of  thermo-mechanical  ero¬ 
sion  a  hypothetic  material,  properties  of  which  are  coincident  with  properties  of 
monofibres.  This  hypothetic  material  is  transversally  isotropic,  as  many  fibres,  in 
particular  carbon  fibres,  are  anisotropic:  properties  along  and  across  fibres  can 
be  very  essentially  different.  Rate  of  thermo- mechanical  erosion  of  fibres  is  deter¬ 
mined  similarly  to  formulae  (72),  (74): 

Def  =  Dfx{p’^,9w){n\  d-na)  +  -^^/||(PE>^tt.)n|.  (73) 

The  concept  of  the  hypothetic  material  of  ’fibres’  allows  us  for  determination 
of  rates  D/j.,  Df\\  to  use  the  models  (developed  for  a  matrix,  therefore  formulae 


for  these  rates  have  the  form  similar  to  (70): 

1/2 


n  =n  ( X(Ps/<^/±)y^^ 


(75) 


where  cr/  and  <j/x  axe  strength  of  monofibres  in  the  longitudinal  and  perpendicular 
directions;  fc/,  p/,  c/  are  heat-physical  characteristics  of  the  fibres;  Jj  and  Eaf 
are  the  pre-exponential  multiplier  and  activation  energy  of  the  process  of  thermo- 
oxidative  decomposition  of  carbon  fibres. 


RATES  OF  EXTERNAL  EROSION  Dj.  and  D\\ 
OF  UNIDIRECTIONAL  COMPOSITES 


Rate  of  external  thermo-mechanical  erosion  D±  in  formulae  (72)  can  be  deter¬ 
mined  with  the  help  of  the  model  presented  in  Fig.  20c  in  the  form  of  a  multilayer 
material  consisting  of  periodically  repeated  layers  of  a  matrix  and  hypothetic  ma¬ 
terial  of  fibres.  Let  thicknesses  of  the  layers  be  hm  aJid  hf,  tm  and  if  be  times 
during  which  there  occurs  a  complete  recession  of  one  layer  of  the  matrix  and  fibre, 
respectively 

tm  =  hmlDue  ,  tf  =  hf/DfA.  (76) 

and  let  us  introduce  ablation  rate  D±  of  the  whole  multilayer  material 


Dj. 


hm  hf 
tm  "1*  if 


(77) 


then  from  (76)  and  (77)  we  obtain: 


D± 


1  ~  hm  ,  hm  \ 

Df±  Djne) 


(78) 


where  hm  =  hml{hf  -|-  hm)  is  the  relative  thickness  of  matrix  interlayer  between 
the  layers  of  fibres  in  the  unidirectional  composite.  It  is  the  desired  formula  that 
connects  ablation  rates  of  matrix  and  fibres  Df±. 
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Rate  of  thermo-mechanical  erosion  on  the  surface  S  being  orthogonal  to  the 
Oxi  axis  under  the  action  of  a  flow  across  the  direction  of  laying  the  fibres  ran 
be  determined  with  the  help  of  the  following  model.  Let  us  consider  a  multilayer 
material  consisting  of  periodically  repeated  layers  of  matrix  and  fibres.  Unlike  the 
case  considered  above,  herein  the  matrix  and  fibres’  layers  are  oriented  across  the 
direction  of  motion  of  the  gas  flow  (Fig.  20d).  Mechanical  erosion  of  matrix  and 
fibres  in  this  multilayer  material  occurs  simultaneously  but  with  different  rates: 
Dm6>  L>/h , therefore  for  the  same  time  </  a  thickness  of  matrix  carried  off  hm 
is  more  than  ft/,  where  hm  =  hf  =  Df\\tf-  There  appears  a  nonsmooth 

surface  S  of  ablation:  ’holes’  forms  in  the  places  of  the  matrix  carried  off.  However, 
when  the  thickness  ft/  reaches  the  certain  value,  recession  of  the  matrix  comes  to 
a  halt,  as  erosional  action  of  the  flow  within  the  ’hole’  eases  towards  its  depth. 
The  recession  process  of  the  matrix  is  resumed,  when  the  thickness  ft/  of  the  fibre 
carried  off  reaches  value  of  hm  (this  model  is  presented  in  Fig.  20c).  Hence,  a 
total  rate  of  recession  of  the  multilayer  material  I?||  will  be  close  numerically  to 
the  value  of  recession  rate  of  the  fibre: 

D\i=Dfn.  (79) 


9.5.4.  FABRIC  COMPOSITES 

Fabric  composites,  for  which  directions  of  warp  and  fill  have  the  same  rights, 
can  be  consider  as  transversally  isotropic  with  a  plane  of  isotropy  being  coincident 
with  the  plane  of  laying  the  fabric.  Let  the  Oxi  axis  be  coincident  with  the  normal 
to  the  plane.  Then  for  such  the  composite,  number  n  =  2  and  tensors  have 
the  following  form  in  the  coordinate  system  Oxi: 

n  =  2,  =  Si28j2  +  SisSjs,  =  SnSji; 

0^2)  _  ^  ^(1)  _ 

and  expressions  (51)  and  (52)  for  rates  of  internal  and  external  erosion  have  the 
form: 


D,  =  f  Dt2{pE,6„)(nl  +nf); 


(81) 
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Di  =  Dii{9^)ni  +  Di2ie^)(nl  +  nl),  (82) 

where  Z)ei  is  the  rate  of  external  erosion  of  the  composite  surface  being  orthogonal 
to  the  Ox\  axis  under  the  action  of  the  flow  onto  the  surface;  De2  is  the  rate  of 
external  erosion  of  the  composite  surface  S  being  orthogonal  to  axis  0x2  or  OX3 
under  the  action  of  the  flow  onto  the  surfaces  respectively,  and  2?||  is  the  rate  of 
internal  erosion  on  the  surface  orthogonal  to  the  Ox  1  axis  and  Di2  -  on  the  surfaces 
orthogonal  to  axes  0x2  and  Oxz . 


RATE  OF  EXTERNAL  EROSION 


To  determine  the  erosion  rate  Dei  let  us  consider  a  model  of  a  multilayer  ma¬ 
terial  consisting  of  periodically  repeated  four  layers  (see  Fig.  20e):  the  flrst  and 
third  layers  are  represented  by  unidirectional  composite  with  perpendicular  (with 
respect  to  the  flow  direction)  laying  its  fibres,  their  recession  rate  is  coincident  with 
D±',  the  fourth  layer  is  an  interlayer  of  the  matrix  between  the  fabric  layers,  its 
erosion  rate  is  Dmc  S’Hd  relative  thickness  is  6'^,  and  the  second  layer  is  considered 
as  unidirectional  composite  turned  by  angle  ‘dm  with  respect  to  the  ags  flow  direc¬ 
tion.  Erosion  rate  of  the  second  layer  D^  is  evaluated  with  the  help  of  formulae 
(72)  and  (73),  in  which  we  should  assume  n2  =  0,  na  =  sini?^,  =  cost9i„: 

=  D||  sin^  ‘dm  +  D±  cos^  dm-  (83) 


A  total  rate  of  recession  of  four- layer  material  Di  is  determined  by  the  method 
of  summing  up,  as  it  was  shown  for  rate  D±  (78)  of  unidirectional  composite: 


Del 


,  Izik  +  iLV'. 

3Dj_  Dm  ) 


(84) 


Erosion  rate  De2  on  fhe  surface  orthogonal  to  the  Oxz  axis  under  the  action  of 
flow  along  the  Oa;i  axis  on  the  basis  of  conclusions  given  above  for  rate  Dy  of 
unidirectional  composite  and  up  to  the  certain  accuracy  can  be  considered  as 
equal  to  rate  Df\\: 


i>e2  =  D||  =  f/ll- 


(86) 


RATE  OF  INTERNAL  EROSION 

Rates  of  internal  erosion  Du  and  Di2  can  be  determined  in  the  similar  way  as 
it  was  made  for  ablative  matrices  with  the  only  difference  that  in  place  of  matrix 
strength  (TmT  we  should  substitute  strength  of  fabric  composite  in  tension  S'^  and 
5^,  respectively: 

(\  ✓  \  2u;— 1 

=  .  (86) 


9.5.5.  BIAS-REINFORCED  FABRIC  COMPOSITES 

Rates  D^i  and  of  internal  and  external  thermo-mechanical  erosion  of 

BRFC  with  angle  V’  in  the  Ox\Xz  plane  can  be  readily  evaluated  with  the  help  of 
formulae  (81)  and  (82)  for  fabric  composite.  In  these  formulae  we  should  assume 
n2  =  0,  ni  =  cos^%  riz  =  sin^; 


De^i  =  Del  cos^  ‘tp  +  De2  sin^  xp;  (87) 

Di^i  =  Dji  cos^  xp  +  Di2  sin^  xp.  (88) 

Total  rate  D^i  of  thermo-mechanical  erosion  of  BRFC  is  determined  as  follows: 

D^i  =  Deijji  -|- Dj^i.+Dc  (89) 


9.5.6.  CALCULATED  RESULTS 

Figure  14  shows  calculated  results  for  different  components  of  thermo-mechanicaI§ 
erosion  of  2-D  composite:  ablation  rate  of  matrix  Dm,  unidirectional  composite 
D||  and  Dj.  and  also  Dei,  De2  and  D^i  -  as  functions  of  surface  temperature 
As  seen  from  the  figure,  the  rate  of  thermo-mechanical  ablation  of  matrix  Dm 
is  essentially  higher  than  the  rate  of  ablation  of  fibres  Dy  and  Dx,  but  rates  of 


ablation  of  fabric  composite  D^i  and  De2  axe  close  to  the  ablation  rate  of  fibres 
in  their  numerical  values. 

Difference  of  rates  Dy  and  D±  is  equal  approximately  to  30  %  under  9y,  =  3000 
K,  and  difference  of  rates  Dei  and  Dt2  of  fabric  composite  in  different  directions 
is  just  higher,  approximately  60  %. 

It  is  explained  by  that  a  matrix  in  fabric  composite  plays  a  greater  role  than  in 
unidirectional  composite. 

As  follows  from  Fig.l4,  magnitude  of  the  internal  erosion  rate  Du  is  also  suffi¬ 
ciently  high:  it  is  equal  approximately  to  25  %  of  the  external  thermo-mechanical 
erosion  rate  Dei- 

Table  13  gives  the  correponding  numerical  values  of  the  functions. 


§9.6.  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
RESULTS  ON  ABLATION  RATE 

To  verify  an  adequacy  of  the  developed  model  of  linear  ablation  we  conducted 
a  comparison  of  numerical  results  of  the  total  linear  rate  of  ablation  Di  on  the 
surface  orthogonal  to  the  Oxi  axis  and  of  experimental  data  for  Generic  2-D 
Caxbon/Phenolic  Fabric  Composites  from  work  [5j.  Velocity  head  was  chosen 
Pa  =  5atm,  pressure  ps  -  equal  to  pj,:  ps  =  Po  =  5  atm.  Rate  Di  was  determined 
in  the  form  of  function  of  surface  temperature 

Di  =  Dc  +  Ds  -l-  Dmi ;  (90) 

Dm\  =  -Del  +  Dii ; 

and  constants  contained  in  the  expression  for  Di 

Jl  E.i,  Ti,  Jj,  E.f,  T,-,  (91) 

5  )  Pm  ?  5  Pf  5  ^0 1 

O'!?  ‘S'l-l-5  •S'2+) 

and  also  Sm  =  0.1,  lo  =  10“^  m 
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were  determined  before  in  [1]. 

It  should  be  noted  that  formulae  (86)  -  (90)  allow  us  to  determine  the  rate  of 
thermo-mechanical  ablation  Dm  of  2-D  composite  with  using  only  mechanical  and 
heat-physical  constants  (91)  without  any  empirical  formulae  for  recession  rates  of 
matrix  and  fibres. 

Figure  15  shows  calculated  results  of  different  rates  of  ablation  for  2-D  car- 
bon/Phenolic  Composite:  external  thermo-mechanical  erosion  rate  D^i,  internal 
erosion  rate  Du  and  also  linear  rate  of  combustion  Dc,  sublimation  Da  and  total 
rate  of  recession  Di.  As  seen  from  this  figure,  under  temperatures  1200  -  2000 
K  a  main  contribution  into  recession  rate  Di  is  of  combustion  rate  Dc  and  exter¬ 
nal  thermo-mechanical  erosion  Z>ei-  Under  temperatures  2000  -  3000  K  diffusive 
combustion  rate  Dc  is  stabilized,  and  Dei  continues  to  grow.  Under  these  temper¬ 
atures  internal  erosion  rate  Du  also  becomes  essential.  Under  temperatures  3000 
-  4000  K  there  occurs  a  sublimation  of  carbon  components  of  the  composite  and 
its  rate  gives  a  main  contribution  to  the  ablation  rate  Di  of  the  composite. 

Points  in  Figure  15  indicate  experimental  values  of  ablation  rate  Di  from  [5]. 
Under  6^^  =  3000  K  experimental  values  of  Di  are  in  the  interval  from  0.5  to 
0.67  mm/s,  and  theoretical  values  are  equal  to  0.6  mm/s.  Thus,  the  accuracy 
of  modelling  of  the  process  of  linear  ablation  of  2-D  Carbon/Phenolic  Composite 
is  sufficiently  high:  it  predicts  the  magnitude  of  Di  falling  within  the  domain  of 
spread  in  experimental  data. 

Figure  16  shows  a  dependence  of  thermo-mechanical  erosion  coefficient  e 
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on  temperature  at  different  values  of  velocity  head  Pq.  It  is  seen,  that  at  initial 
period  of  heating  under  900  -  1200  K  values  of  e  are  the  closest  to  1,  as  recession 
is  determined,  in  the  main,  by  mechanical  ablation.  Under  temperatures  1200  - 
2000  K,  as  mentioned  above,  an  essential  contribution  is  of  combustion  rate  D^ 
therefore,  values  of  e  decrease  to  0.4  -  0.6.  Since  under  temperatures  2000  -  3000 
K  diffusive  combustion  rate  Dc  is  stabilized,  value  of  e  anew  grows  up  to  0.6  - 
0.8.  And,  finally,  when  sublimation  process  starts  within  the  interval  9^  =  3000 
-  4000  K,  value  of  the  thermo-mechanical  erosion  coefficient  falls  down  practically 
to  zero. 
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Fig.  14. 

Varying  Rates  of  Thermo-mechanical  Recession  of 
Phenolic  Matrix  (Dm),  Unidirectional  Composite  {Dj)  and  (T)//)  and 
2-D  Fabric  Composite  {Dei,De2)  and  Internal  Erosion  Rate  Da 
of  2-D  Carbon/Phenolic  Composite  versus  temperature  9  at  p^=  5.25  atm. 
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Rates  of  Thermo -Mecanical  Ablation  : 

1. Exteranl  T-M  Erosion  for  Phenolic  Matrix  , 

2. Exteranl  T-M  Erosion  for  Unidirectional 
Carbon/Phenolic  Composite  along  Fibers  Du  , 

S.Exteranl  T-M  Erosion  for  Unidirectional 

Carbon/Phenolic  Composite  Across  Fibers  Dj^  , 

4. Exteranl  T-M  Erosion  for  2-D  Textile 

Carbon/Phenolic  Composite  along  Fabric  , 

5. Exteranl  T-M  Erosion  for  2-D  Textile 

Carbon/Phenolic  Composite  Across  Ply  Fabric  D^^  , 

6. Internal  T-M  Erosion  for  2-D  Textile 

2-D  Textile  Carbon/ Phenolic  Composite  Oh  . 


Temperature,  K  Rates  of  Recession,  mm/s 
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Fig.  15. 

Varying  Rates  of  Thermo- mechanical  {Dei  and  and 
Thermo-chemical  Ablation  {Dc  and  Z?,)  and  total  rate  of  ablation  Di 
of  2-D  Carbon/Phenolic  Composite  versus  temperature  0 
at  p  =  5.25  atm  in  air.  Solid  curves  -  calculated  results, 
points  -  experimental  data  of  JDi  [5^]. 
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Relation  Between  Rates  of  Thermal -Mechanical  and 
Thermo -chemical  Ablation  of  2-D  Carbon/ Phenolic 
Composite  at  P=5  atm  in  air 


Temperature  Rate  of  Recession ,  mm/s 


Pur  K 


Subli - 

mat ion 

Combus - 

tion 

Dc 

External 

T-M  erosion 
^ei 

Internal 

T-M  erosion 
^ii 

Total 

1500.0 

.000 

.047 

.077 

.018 

.142 

2000. 0 

.000 

.127 

.154 

.049 

.329 

2500.0 

.000 

.135 
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.002 
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.196 
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4000. 0 
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.254 

2.331 
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Fig.  16. 

Varying  Coefficient  of  Thermo-mechanical  Erosion  € 
of  2-D  Carbon/Phenolic  Composite  versus  temperature 
at  different  values  of  valocity  head  in  air. 
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Varying  Coefficient  of  Thermal -Mechanical 
Erosion  of  2-D  Carbon/Phenolic 
Composite  on  Temperature  at  different 
values  of  Velosity  Head  in  air 
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Fig.  17. 


Varying  Ablation  Rates  of  Carbon/Phenolic  Composite 
versus  temperature  9  at  different  bias  angles  xj) 
and  Po  =  3  atm  in  air. 
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Varying  of  Total  Ablation  Rates  of  Carbon/Phenolic  BRCM 
vs  Temperature  for  different  values  of  angle  ^  at  P  =3  atm 
in  air 


Temperature f  K 

e,.. 

Rate  of 

Recession 

mm/s 

angle  =0® 

T^30'^ 
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y^=6(f 

1500.0 
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.531 
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.671 
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.491 
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With  increasing  the  velocity  head  Po,  value  of  the  coefficient  e  also  grows  (see 
Fig.16). 

The  corresponding  numerical  values  are  given  in  Table  15. 

Figure  17  shows  calculated  results  of  the  total  rate  of  recession  D^i  of  composite 
BRFC  for  Po  =  3  atm  (ps  =  pj,  =  3  atm)  as  a  function  of  surface  temperature  6^ 
for  different  values  of  angle  ?/>  =  0,  15,  30,  45,  60,  75°. 

Heat-transfer  coefficient  (a/cp)  was  assumed  to  be  equal  to  (a/cp)  =  1.1  kg/m^- 
s.  As  seen  from  Figure  17,  with  increasing  angle  '0  the  recession  rate  of  BRFC 
lowers  within  all  the  temperature  interval.  Under  =  3000  K  values  of  rate 
D^i  are  equal  to  0.53,  0.48  and  0.34  for  angles  0°,  30°  and  75°,  respectively.  Thus, 
for  small  angles  0  (0°  -  30°)  lowering  the  ablation  rate  is  not  very  essential  («  10 
%  as  compared  with  angle  0  =  0°),  and  only  for  greater  angles  (60°  -  75°)  this 
lowering  is  equal  approximately  to  35  %. 

The  corresponding  numerical  values  D^i  are  given  in  Table  6. 


§9.7.  CALCULATION  OF  HEAT  FLUX 
ON  ABLATIVE  SURFACE 

To  determine  true  temperature  of  ablative  surface  6^  of  composite  under  the 
action  of  overrunning  hot  gas  flow,  it  is  necessary  to  take  account  of  processes  of 
heat  and  mass  transfer  occuring  on  this  surface. 

Total  heat  flux  to  the  ablative  surface  of  the  composite,  which  is  consumed 
for  warming  up  the  material,  consists  of  the  following  components: 

9e  —  9e  ~  QR  ~  ~  Qw  ~  Qsub  “I"  9cs  "H  Qcgi  (^3) 


where 

=  [afcp)(Hr  —  Ho)  -  convective  heat  flux  to  the  ’cold  wall’  supplied  from  the 
heated  boundary  layer  to  the  composite  surface  S, 

=  (afcp)(Hjji,  —  Ho)  -  convective  heat  flux  carried  off  from  the  surface  S  heated 
up  to  temperature  into  the  surroundings, 

Qfi  =  -  radiational  heat  flux  carried  off  from  the  surface  S  at  the  expense 

of  own  radiation  of  the  surface. 


Qbi  =  lbi{pg'^g)w{Hr  —  Hyj)  -  convective  heat  flux  carried  off  from  the  heated 
surface  S  at  the  expense  of  injection  into  the  boundary  layer  of  gaseous  products 
of  thermo-decomposition  of  phenolic  matrix, 

<lsub  =  PwDgAQaub  ’  heat  flux  carried  off  from  the  composite  surface  S  at  the 
expense  of  physico-chemical  process  of  sublimation  of  the  composite, 

Qcs  —  PwDcY\AQcOjYI  -  heat  flux  supplied  to  the  composite  surface  S  at  the  ex¬ 
pense  of  exothermal  process  of  combustion  of  solid  phases  of  the  composite  in  the 
oxidizing  overrvmning  flow, 

<lcg  =  'ybi(PgVg)wDc  AYiAQi  -  heat  flux  supplied  to  the  composite  surface  S  at 
the  expense  of  exothermal  process  of  combustion  of  gaseous  products  of  thermo¬ 
decomposition  of  polymer  binder  injected  into  the  boundary  layer. 

There  is  a  following  notation: 

(a/cp)  -  heat-transfer  coefficient  between  the  composite  surface  E  and  boundary 
layer  of  hot  gas; 

Hr  -  enthalpy  of  the  overrunning  flow; 

Hup  =  cdw  -  enthelpy  of  the  composite  surface  heated; 

Ho  =  c6o  -  enthalpy  of  the  ’cold  wall’; 
asB  -  Stefan-Boltsman  constant; 

Ew  -  degree  of  blackness  of  the  heated  surface; 

76/  -  injection  coefficient  (0  <  75/  <  1); 

AQsub  -  heat  of  sublimation  of  carbon  phases  of  the  composite; 

AQ 

com  "  heat  of  combustion  of  solid  phases  of  the  composite  (absolute  value); 

AQi  -  combustion  heat  of  gaseous  products  of  phenolic  matrix  thermo-decomposition| 
(absolute  value); 

Y\  -  mass  concentration  of  carbon  component  on  the  composite  surface  S; 

AYi  -  changing  mass  concentrations  of  chemical  components  in  the  mixture  of 
gaseous  products  of  thermo-decomposition  of  phenolic  matrix  at  combustion  on 
the  surface  S; 

{PgVg)w  =  —{KVpg)w  •  n  -  mass  rate  of  injection  of  gaseous  products  of  phenolic 
matrix  thermo-decomposition  into  the  boundary  layer, 

C^Pg)w  -  pore  pressure  gradienton  the  ablation  surface  E. 


To  conduct  numerical  calculations,  we  use  the  following  values  of  constants: 


(a/cp)  =  1.1  kg/m^  •  s;  =  12  MW/m^;  6o  =  296K; 

£xd  —  0.8;  ~  0.4;  ^Qsub  ~  20  MJ/kg; 

AQeom  =  9.8  MJ/kg;  Yi  =  1. 


Fig.  18. 

Dependence  of  equilibrium  temperature  Or  and 
true  temperature  of  the  surface  6^  (•) 
of  2-D  Carbon/Phenolic  Composite  on  Cold  Wall  Convective  Heat  Flux  qt. 
Solid  curves  -  calculated  results,  points  -  experimental  datat  5  J . 
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Varying  of  Equilibrium  and  True 
Temperatures  of  Ablation  Surface 
vs  Cold  Wall  Heat  Flux  ^ 


Heat  flux,  MW/m2 


Temperature 
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Temperature  ^  Temperature 
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Fig.  19. 

Varying  Rates  of  Thermo-mechanical  {D^,i  -t-  Di^, 
Thermo-chemical  {Djc)  and  total  Di  Ablation 
of  2-D  Carbon/Phenolic  Composite  versus  Heat  Flux  at  pj=  5  atm 
in  air.  Solid  lines  -  calculated  results,  •  -  experiments  [5’]. 
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Varying  RatBS  of  Th&naal-MBchanical  and 
Thermo -Chemical  Ablation  of  2-D  Carbon/Phenolic 
Composite  on  Heat  Flux  at  P=5  atm  in  air 


HEAT  FLUX,  MVf/m2  Rate  of  Recession,  mm/s 
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Chapter  X 

CALCULATION  OF  INTERNAL  HEAT-MASS-TRANSFER 
AND  THERMOSTRESSES  IN  A  CYLINDRICAL  SHELL 
MADE  OF  BRFC  UNDER  THE  ACTION 
OF  A  HIGH-TEMPERATURE  GAS  FLOW 

The  aim  of  the  present  Chapter  is  to  conduct  calculations  of  fields  of  tem¬ 
perature  9,  pore  gas  pressuure  p,  volumetrical  phase  concentrations  cpj,  and  also 
thermostresses  caused  by  nonuniform  heating  and  volumetrical  ablation  in  a  cylin¬ 
drical  shell  (Figure  1)  made  of  Bias- Reinforced  Carbon/Phenolic  Composite  under 
the  action  of  a  high-temperature  gas  flow  onto  the  shell. 


§10.1.  MATHEMATICAL  PROBLEM  STATEMENT 

A  general  statement  of  the  mathematical  problem  on  coupled  calculation  of 
internal  heat-mass-transfer  and  thermostresses  in  ablative  composite  materials  has 
been  formulated  in  [1]  and  also  in  our  previous  works  [6-9].  Let  us  apply  this  model 
for  calculation  of  a  cylindrical  shell  made  of  BRFC  (Figure  1). 

Let  us  consider  further  a  cylindrical  coordinate  system  Ordz  and  assume  that 
the  axes  of  the  system  are  coincident  with  the  axes  of  the  coordinate  system  Oxi 
considered  above  and  rotated  by  angle  V’  about  principal  orthotropy  axes  Ox[  of 
the  composite. 

Let  us  consider  a  hollow  cylinder  of  ablative  material  BRFC  subjected  to 
intensive  heating  on  the  external  side.  Convective  heat  fluxes  and  pressures  pe 
are  uniformly  distributed  on  the  lateral  surface. 

There  is  a  total  pressure  of  overrunning  flow  pe  =  peo  +  Po  at  the  external 
surface  of  the  cylinder  r  =  Rq,  where  peo  is  the  static  pressure  of  quiescent  at¬ 
mosphere,  Po  is  the  velocity  head  of  gas  p[,  =  Pgevl^,  Pge  is  the  density  of  gas 
flow  on  the  composite  surface,  =  V3  is  the  component  of  gas  velocity  V  at  the 
external  surface  of  the  boundary  layer  (in  the  Oz  axis  direction).  Projection  of 
the  flow  velocity  vector  V  on  the  cylinder  surface  r  =  Rd  (ablation  surface  E)  is 


-  BO- 


determined  from  (46)  as:  Vei  =  Vs2  =  0;  V^z  =Vz=  v^z-  Pressure  ps  from  (50) 
is  determined  as  follows:  ps  =  PgeVz  —  Po-  internail  surface  of  the  cylinder 

r  =  Ri  a  pressure  is  assumed  to  be  equal  to  static  pressure  peo- 

Both  the  ends  of  the  cylindrical  shell  z  =  0  and  z  =  L  are  assumed  to  be  free 
from  stresses,  i.e.  cTz  =  — Peo  at  z  =  0,  L.  The  ends  are  assumed  to  be  hermetic 
and  heat- insulated. 

Let  us  state  the  following  problem:  to  determine  parameters  of  internal  heat- 
mass-transfer  in  an  ablating  structure  and  to  find  stresses  in  the  structme  caused 
by  external  heating  and  ablation  (volumetric  and  surface). 

Conditions  of  the  problem  allow  to  use  a  statement  in  terms  of  stresses  (see 
item  [1]). 

In  the  scope  of  these  assumptions  a  solution  of  the  problem  is  sought  in  the 
form  of  functions  depending  on  radius  r  and  t.  Stresses  are  determined  as  follows: 

F 

(Tr  = - h  <PoP;  cr^  F'  +  <Pgp-,  (Trz  =  0;  (94) 

r 


—  PeOf 

where  F  is  the  stress  function,  and  F'  =  dF/ dr. 

Constitutive  relations  (12)  from  [1]  have  the  following  form  for  the  considered 
problem: 


O  crr  +  flil,p  1^X2^  (  ,  7  \  ViZtjj  .  7  . 

£r  —  Sllii  H - - S  +  J2rt>p)  s  W3V>p  —  PeOJi 

JiilTP  Jhl 


o  Vi2^  j-  .  CTj)  -f /2^p  V2Z  f7  \  (c.r.\ 

=  S22ip  —  +  fltl^p)  H - ^  -^r~\h^P  ~  PeOjt 


2i) 


E- 


2^ 


E2ij} 


+  flipp)  ~  jp  ^(<^1?  +  f2xpp)  +  „  {fzipp  Peo) 


^'23V> , 


E 


Ixjj 


E, 


2xp 


E: 


3^ 


Elasticity  moduli  Ei^  and  Poisson  ratios  Uij^  are  expressed  by  formulae  (10), 
heat  deformations  entp  -  by  formulae  (7)  and  coefficients  of  phase  interaction  fi^ 
-  by  the  formulae: 


firp  =  fi  cos^  ^4-/3  sin^  V*; 
fz^  =  /i  sin^  V’  +  /a  cos^ 

/2^  =  /25 


(96) 


where  functions  fi  depend  on  composite  porosity  ipg  and  are  expressed  as  follows: 


fi  =  fm  = 


-  <p!) 


f2  —  ^C(l  “  ^/^)  +  1  “  /m! 

/3=((i-c)(i-yrf)+c)/™. 


(97) 


Here  (f°  is  the  volumetric  concentration  of  fibres  in  the  composite,  (  is  the 
ratio  of  the  fibres’  number,  per  unit  volume  of  the  composite,  directed  along  the 
weft  to  the  fibres’  number  directed  along  the  fabric  warp. 

To  determine  function  F  we  have  the  equation  of  compatibility  of  deforma¬ 
tions: 

de^  —  Er 


+ 


=  0, 


dr  r 

into  which  we  should  substitute  relation  (95).  As  a  result,  we  obtain  the  equation 
of  the  second  order  for  F\ 


d^F 

dr'^ 


d 


+  vr- 


dr  \E. 


9F 

dr 


rE 


Irjj 


F\^r 


+  I  1  £2#  +  E2*|:  ((-<^12*7i*  +  In,  -  ■'23*7*)  ^)  + 


dr 


+  -  (-I^12V'/lV>  +  h^l)  -  V2zh%l>  ~  ifl^  -  i'l2V’/2^  -  =  0. 

r  \  j 


Heat-mass-transfer  equations  have  the  form: 

d(pa 


(98) 


Pi 


dt 


n  T 


d  \  d  f  d  \  ^ 

-^{pg'pg)  =  {^^^ii^'^i.^gPg^))  +  -1- J/F/; 

^  ^  1=1 


(99) 
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Boundaxy  and  initial  conditions  for  the  set  of  equations  (98),  (99)  are: 


r  =  Ri  : 

F  =  —RlPel\  RpgO^Pel\ 

(100) 

r  =  Rd{J)  '■ 

F  =  RDpe2’>  Rapg^  ~  Pe2?  ^11^  “  9e2  5 

<  =  0: 

<Pa  “  ^1?  Vi  “  yiO?  Pg  ^  Pg^^  9  — 

(101) 

here  Ri  is  the  internal  radius  of  the  cylinder,  Roif)  is  the  external  radius  decreas¬ 
ing  with  time  due  to  linear  ablation  of  the  composite: 

RD  =  R2-f  D^i{r)dT,  (102) 

Jo 

where  R2  is  the  external  radius  of  the  cylinder  at  initial  time  t  =  0,  is  the 
total  rate  of  linear  ablation  of  BRFC  calculated  by  formula  (89). 

Heat  flux  5ei  carried  off  from  the  internal  surface  r  =  Ri  of  the  shell  is 
determined  by  convective  heat  transfer  and  ragiation: 

qel=  OCTl{0  — 6 el)  SB (103) 

here  an  is  the  coefficient  of  heat  transfer  from  the  internal  surface,  9ei  is  the 
temperature  of  the  surrounding  air  inside  the  cylinder,  asB  is  the  Stefan-Boltsman 
coefficient,  e  is  the  degree  of  blackness  of  the  internal  surface. 

Heat  flux  5e2  supplied  to  the  external  recessive  surface  r  =  Roi^)  of 
composite  shell  is  calculated  by  the  relations  (93): 

qe2  =  9e-  (104) 

The  equation  system  (98)  -  (99)  with  conditions  (100)  -  (101)  is  the  mathe¬ 
matical  statement  of  the  problem  on  determination  of  functions  F,  (p^,  Vi,  Pg  and 
6  depending  on  coordinate  r  and  time  t.  On  determining  a  value  of  stress  function 
jP,  stresses  (Tr  and  in  the  cylindrical  composite  shell  can  be  found  by  formulae 
(94). 

Pore  gas  pressure  p  is  determined  by  the  formula: 


P  —  RaPg^' 


(105) 
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§10.2.  DAMAGE  PARAMETERS 


In  accordance  with  the  model  developed,  fabric  ablative  composites  and  BRFC| 
have  four  main  types  of  destruction  (A),  (B),  (C)  and  (D)  (see  [1]).  For  the  cylin¬ 
drical  shell  made  of  BRFC  two  types  of  destruction  are  possible:  (A)  -  breakage 
of  threads  laid  in  Warp  and  Fill  Directions  of  Fabric  0x\  and  Ox'^  and  (C)  - 
delamination  of  the  composite. 

Destruction  of  the  type  (A)  starts  when  damage  parameter  Z2{r,t)  reaches 
value  1  within  a  certain  zone  with  coordinate  r  =  at  time  t^: 


Z2{riq)  =  l- 


(106) 


this  parameter,  from  relations  (14)  of  [1]  and  equations  (14),  (16)  of  the  present 
work,  is  determined  by  the  formula: 


c+ 


(107) 


As  it  was  noted  in  Chapter  VIII,  for  angles  V’  >  5°  the  condition  <  1 

is  usually  satisfied,  therefore  from  (107)  it  follows  that 


Z2{r,t)  =  — qj-  at  Ip  >  5°. 


(108) 


At  the  same  time  for  fabric  2-D  composites  with  angle  tp  =  0°  the  condition 
>  1  is  satisfied,  and,  thus,  damage  parameter  Z2  takes  the  form: 


.  .  (7^  cos^  ‘tp 

Z2{r,t)  =  — -  at  0  =  0 


(109) 


These  expressions  for  Z2  at  different  angles  xp  will  be  used  in  calculations  below. 

Destruction  by  the  type  (C)  starts  when  damage  parameter  zi{r,t)  reaches 
value  1  within  a  certain  zone  of  the  cylinder  with  coordinate  r  =  at  time 


=  (110) 

which  from  equations  (15)  of  [1]  and  (12)  of  §8.2  is  determined  by  the  formula: 

(Trir,t) 


zi{r,t)  = 


(111) 


Expression  (111)  is  valid  for  all  angles  V’  from  0°  to  90°,  and  at  ^  =  0°  we  have 
5+  -  5+ 


$  10.3.  METHOD  OF  NUMERICAL  SOLUTION 
OF  THE  PROBLEM 

The  peculiarity  of  the  problem  statement  given  above  consists  in  that  a  do¬ 
main  of  its  resolution  is  variable:  it  is  interval  Ri  <  r  <  Ro{t)-  For  numerical 
solving  the  problem  within  variable  domains  we  need  to  reconstruct  difference 
schemes  or  finite-element  grids  at  each  time  step  that  usually  leads  to  an  appear¬ 
ance  of  calculation  instability  and  lowering  the  accuracy  of  calculated  results.  To 
avoid  this,  we  developed  the  method  of  passage  to  mobile  coordinates  [6].  For  the 
problem  considered,  this  method  consists  in  the  following. 

Let  us  introduce  a  mobile  coordinate  x: 


X  —  a{t)r  -1-  6(t), 

where  functions  a(t)  and  b{t)  are  determined  by  the  formulae: 

_  R2  -  Ri  _  RiiR2  —  Rpii)) 

'  ~  Ro(t)  -Ri'  ^  Roit)  - Ri 


(112) 


(113) 


The  transformation  of  coordinates  (112),  (113)  reduces  the  mobile  domain 
Ri  <  r  <  Roit)  to  the  non-mobile  Ri  <  x  <  R2- 

The  equation  system  (98)  -  (99)  in  the  mobile  coordinate  system  is  written 
as  follows: 

‘^a  ~  ~  f I Pf  ’i 

Vi  —  —  —JilP2'i 

2  ^ 
i^gPgY  -  Drl’l^iV’gPg)'  =  ^  {RaKii{x  -  h){pg9)')  +  ^ 


pcd  -  pcD^^e'  = - rihi^ix  -  h)0')'  +  Cga^R  -  aIU-^9' {pg6)' - 

X  —  0 


VpSAe”£j,ri-Ae5v>!J/r/, 


(114) 


and  compatibility  equation  is: 


X  —  b  El 


Xx  - 


E2iI}F  -\- 


+  1  + 


£22t/>  ~ 


p!'  ,  ( _1 _ F'  -i  — _ -  f  1  ^  E2i,F-\- 

Vx-6  £^2V’/  y(x-6)Ei^  \(x-6)£?i^/y 

+  (^22ip  +  ^  ^-^'12^/lV’  +  ~  ^23rpf3ip^^  + 

H - ^  ^^'l2v/lV'  d"  f2ip  ~  ^23ipf34>  ~  ^  ~  ^12il)f2il)  ~  =  0- 

Here  derivative  with  respect  to  coordinate  x  is  (')  =  d/dx.  Boundary  (100)  and 
initial  (101)  conditions  remains  their  forms  in  coordinates  x,t  excepting  the  heat 
conditions: 

X  —  Rx  :  cik\xip^  —  Qci ^ 


X  =  Ro 


Clkxixl^O  —  ^e2  • 


(116) 


Expressions  for  stresses  are  transformed  in  the  similar  way; 

aF 

CTr  =  +  ‘fgP, 

=  aF'  +  (fgp.  (117) 

Additional  function  ^  =  ^(x,t)  is  introduced  by  equation: 

^  R2  -  Rx  { _p\  (llg) 

^  {RD-Rif\  a 

Here  all  unknown  functions  F^  /?p,  <^05  Vi  depend  on  arguments  x  and  t. 

Equation  system  (114)  -  (116)  has  been  solved  numerically  with  the  help  of 
step-by-step  method  and  difference  implicit  absolutely  stable  scheme.  To  solve  the 
systems  of  linear  algebraic  equations  arising  in  applying  the  difference  method,  we 
use  the  procedure  of  matrix  sweeping. 


(118) 


§10.4.  CALCULATED  RESULTS. 
INITIAL  DATA. 


In  accordance  with  the  method  described  above,  numerical  calculations  were 
conducted  to  determine  a  thermomechanical  behaviour  of  a  cylindrical  shell  made 


of  Carbon/phenolic  BRFC.  Characteristics  of  the  composite  were  given  the  same  as 
they  were  presented  in  chapter  VII  -  IX  of  the  present  report  and  in  [1],  Conditions 
of  heating  were  chosen  according  to  [11]: 

heat  flux  to  the  cold  wall  5°  =  12  MW/m^; 

total  pressure  of  gas  flow 

on  the  external  surface  of  the  shell  Pe2  =  3  atm; 

pressure  on  the  internal 

surface  of  the  shell  Pei  =  1  atm. 

Besides  that,  we  assume  the  following  values  of  heat  characteristics  of  the 

flow: 

heat-transfer  coefficient  with 
hot  gas  flow 
heating  time 

chemical  composition  of  gas 
at  external  surface  of 
boundary  layer 
heat-transfer  coefficient  from 
internal  surface  of  the  shell 
degree  of  blackness 

Geometrical  dimensions  of  the  shell: 

thickness  12  mm; 

internal  radius  125  mm; 

Bias  angle  tp  0°  and  30°. 

Numerical  calculations  were  conducted  for  two  angles  •tp  =  0°  and  30°. 


{oi/cp)  —  2.2  kg/m^  •  s; 
t  =  7...10  s; 
corresponds  to  standard 
composition  of  air  at 
d  =  295A' 

oi'^  =  10  W/(m2  .  K); 

£u)  =  0.8. 


10.4.1.  THE  EFFECT  OF  ABLATION  ON  THERMOMECHANICAL 
BEHAVIOUR  OF  A  SHELL  OF  BRFC 

At  first  let  us  describe  results  for  V’  =  30° . 

Figure  21  shows  changing  characteristics  of  BRFC  corresponding  to  these 
conditions  of  heating:  temperature  9w  of  ablating  surface  r  =  RDi't)^  linear  rate 
of  ablation  D^i  and  temperature  9a  of  the  internal  surface  of  the  shell  r  =  Ri 
depending  on  time  t  of  heating.  As  seen  from  Fig.21,  temperature  9xu  of  the 
ablating  surface  quickly  (during  0.1  -  0.2  s)  comes  to  a  quasistationary  regime  and 
does  not  practically  change  during  all  the  time  of  heating  being  at  the  level  of 
3000  K.  The  same  character  has  function  D^i(t):  its  quasistationary  magnitude 
is  equal  to  0.62  mm/s  for  angle  ^  =  30°  under  the  given  conditions  of  heating. 

Temperature  of  the  internal  surface  9s  starts  to  differ  essentially  from  initial 
temperature  =  295  K  only  in  f  6  s  and  at  time  t  =  7  s  is  equal  to  414  K 
(Fig.21). 

Figure  22  shows  distributions  of  temperatme  9(r,  t)  versus  the  shell  thickness 
for  different  times  t  =  0.5,  1,  3,  5  and  7  s.  As  seen  from  this  figure,  during  heating 
time  t  the  external  surface  of  the  shell  r  =  Roit)  is  carried  off  and  temperature 
profile  9[r,  t)  is  displaced  in  the  direction  of  internal  surface  r  =  Ri . 

Corresponding  distributions  of  pore  gas  pressure  p{r,t)  versus  the  shell  thick¬ 
ness  are  given  in  Figure  23.  As  seen  from  this  figure,  the  level  of  pore  gas  pressure 
for  the  shell  with  angle  ^  =  30°  is  sufficiently  high  and  is  equal  to  Pmax  =  40  — 100 
atm.  Distribution  of  pore  pressure  p(r,  t)  versus  coordinate  r  has  a  typical  peak, 
maximum  of  which  Pmax  corresponds  to  the  temperature  value  9{r^i)  of  the  shell 
and  9{r,  t)  w  600  K  (see  Figs.  22  and  23).  This  peakis  displaced,  as  the  composite 
is  warmed  up,  towards  the  internal  surface  r  =  Ri  of  the  shell,  and  a  maximiun 
value  of  pore  pressure  Pmaxif)  increases  with  time  from  40  atm  at  f  =  0.5  s  to  ~ 
100  atm  at  <  =  7  s. 

Distributions  of  volumetric  concentration  of  polymer  phase  of  the  matrix 
<P2{r,t)  versus  the  shell  thickness  are  shown  in  Figure  24  for  different  times.  As 
seen  fron  thi  figure,  within  the  warmed  zone,  where  9{r,  t)  >  600  K,  value  of 
sharply  falls  down  from  initial  value  =  0-35  to  0.  Value  (p2{r,t)  =  0  corresponds 
to  the  domain  of  composite  coked  completely.  The  width  of  the  domain  during 


-ss- 


(0,^00  >  '0'V 


Dttv  mm/s 


time  t,  s 


Fig.  21. 

Varying  external  surface  temperature  9^, 
temperature  Os  on  the  internal  surface  R  =  r\  and 
recession  rate  i^js^ersus  duration  of  heating  t  of  the  shell 
made  of  BRFC  with  0  =  30°. 
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Tabl.  19 

Varying  of  True  Temperature  of  Surface , 
Rate  of  Recession  and  Temperature  of 
Internal  Surface  of  BRFC  vitliy''=30® 
vs  Time  of  Heating 


Time  t,s  Temperature 

^ 


Rate  of  Temperature 

Recession  P$  ,K. 

Dtfrj  mm/s 


.1 

2143.0 

.2 

2911.6 

.3 

2972.7 

.4 

2988.9 

.5 

3004.3 

.6 

3014.1 

.  7 

3024.6 

.8 

3030.9 

,9 

3036.5 

1.0 

3041.3 

1.2 

3049.1 

1.4 

3055.4 

1.6 

3060.5 

1.8 

3066.5 

2.0 

3070.2 

2.2 

3073.2 

2.4 

3076.0 

2.  6 

3078.3 

2.8 

3080.5 

3.0 

3082.4 

3.2 

3084.2 

3.4 

3085.8 

3.6 

3087.3 

3.8 

3088.6 

4.0 

3092.1 

4.2 

3092.9 

4.4 

3094.0 

4.6 

3095.0 

4.8 

3095.9 

5.0 

3096.8 

5.2 

3097.6 

5.4 

3098.4 

5.6 

3099.1 

5.8 

3099.8 

6.0 

3100.5 

6.2 

3101.1 

6.4 

3101.7 

6.6 

3102.3 

6.8 

3102.9 

7.0 

3103.4 

7.2 

3103.9 

7.4 

3104.4 

.000 

294.0 

.390 

294.0 

.602 

294.0 

.620 

294.0 

.625 

294.0 

.630 

294.0 

.633 

294.0 

.636 

294.0 

.638 

294.0 

.640 

294.0 

.643 

294.0 

.645 

294.0 

.647 

294.0 

.649 

294.0 

.650 

294.0 

.652 

294.0 

.652 

294.0 

.653 

294.0 

.654 

294.0 

.655 

294.0 

.655 

294.0 

.656 

294.0 

.657 

294.1 

.657 

294.2 

.657 

294.4 

.658 

294.7 

.659 

295.3 

.659 

296.2 

.660 

297.8 

.660 

300.3 

.660 

304.1 

.660 

309.  7 

.661 

317.7 

.661 

329.1 

.661 

344.6 

.661 

365.5 

.662 

393.0 

.662 

428.6 

.662 

473.  7 

.662 

530.0 

.662 

599.4 

.663 

684.2 

~30 


Tahl.  20 

Varying  of  True  Temperature  of  Surface ^ 
Rate  of  Recession  and  Temperature  of 
Internal  Surface  of  BRFC  with  ^^=0® 
vs  Time  of  Heating 


Time  t,s 

Temperature 
,  K 

Rate  of 

Recession 

D^jmm/s 

Temperature 

Og  ,R 

.1 

2143.  0 

.000 

294.0 

.2 

2911.6 

.363 

294.0 

.3 

2972. 7 

.558 

294.0 

.4 

2988.9 

.574 

294.0 

.5 

3004.3 

.579 

294.0 

.6 

3014.1 

.583 

294.0 

.  7 

3024. 6 

.586 

294.0 

.8 

3030.9 

.589 

294.0 

.9 

3036.5 

.591 

294.0 

1.0 

3041.3 

.593 

294.0 

1.2 

3049.1 

.595 

294.0 

1.4 

3055.4 

.598 

294.0 

1.6 

3060.5 

.599 

294.0 

1.8 

3066.5 

.601 

294.0 

2.0 

3070.2 

.603 

294.0 

2.2 

3073.2 

.604 

294.0 

2.4 

3076.  0 

.604 

294.0 

2.  6 

3078.3 

.605 

294.0 

2.8 

3080.5 

.606 

294.0 

3.0 

3082.4 

.607 

294.0 

3.2 

3084.2 

.607 

294.0 

3.4 

3085.8 

.608 

294.0 

3.6 

3087.3 

.608 

294.0 

3.8 

3088. 6 

.609 

294.0 

4.0 

3092.1 

.609 

294.0 

4.2 

3092.9 

.610 

294.0 

4.4 

3094.0 

.610 

294.0 

4.6 

3095. 0 

.611 

294.0 

4.8 

3095.9 

.611 

294.0 

5.0 

3096.  8 

.611 

294.0 

5.2 

3097.  6 

.612 

294.0 

5.4 

3098.4 

.612 

294.0 

5.6 

3099.5 

.612 

294.0 

5.8 

3099. 8 

.612 

294.0 

6.  0 

3100.5 

.613 

294.1 

6.2 

3101.1 

.613 

294.2 

6.4 

3101.7 

.613 

294.3 

6.  6 

3102.3 

.613 

294.6 

6.8 

3102.9 

*  .613 

295.1 

7.0 

3103.4 

.614 

295.9 

7.2 

3103.9 

.614 

297.4 

7.4 

3104.4 

.614 

299.8 

T—  R  j  P—  R2 


Fig.  22. 


Distribution  of  temperature  d  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t  =  0.5,  1,  3,  5,  7  s 
(r  -  radius  of  the  shell).  Symbol  <—  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Roit)  at  times  t.‘ 


1.25  1.27  1.29  1.31  1.33  1.35  r 


r=Ri 


r=R2 


Fig.  23. 


Distribution  of  pore  gas  pressure  p  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t 
(r  -  radius  of  the  shell).  Symbol  <—  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Roit)  at  times  t. 


warming  up  grows  and  becomes  equal  to  3.5  mm  at  t  =  7  s,  while  the  thickness  of 
the  composite  layer  carried  off  is  equal  to  4  mm  at  the  same  time  moment. 

Figure  25  shows  distributions  of  tangential  stresses  f)  versus  the  shell 
thickness  at  different  times  t  of  heating.  At  initial  times,  when  t  <  0.2  s  and  ab¬ 
lation  process  of  the  composite  (volumetric  and  surface)  is  not  yet  realized,  heat 
deformation  £22^  of  the  composite  is  positive,  that  leads  to  the  appearance  of 
compressive  (negative)  values  of  stresses  cr^.  As  the  composite  is  warmed  up,  for 
t  >  0.5  s  near  the  external  surface  r  —  Rd  sl  volumetrical  ablation  starts  that  leads 
to  the  appearance  of  negative  values  of  heat  deformations  £22^/'!  thus  there  appears 
a  shrinkage  which,  in  its  turn,  causes  an  appearance  of  tensile  (positive)  values  of 
tangential  stresses  .  For  the  case,  within  the  domain  of  relatively  low  tempera¬ 
tures,  where  9(r,t)  <  1500  K,  there  are  compressive  stresses  This  stress  profile 
<7^{r,t)  is  displaced  towards  the  internal  surface  of  the  cylinder,  as  the  composite 
is  warmed  up.  A  maximum  level  of  tensile  stresses  grows  max  60  MPa. 
However,  further,  near  the  external  ablating  surface,  stresses  start  to  decrease 
(approximately  from  <  =  2  s)  and  go  anew  into  the  domain  of  negative  values  (at 
t  —  7  s  a^{Ro,t)  —  —15  MPa).  This  phenomenon  is  caused  by  peculiarities  of 
heat  deformation  £22^^  of  ^he  type  of  Carbon/Phenolic  composite,  for  which  under 
temperatures  6  >  2400  K  the  shrinkage  is  replaced  by  heat  expansion  that  leads 
to  the  appearance  anew  of  negative  values  of  within  the  domain  adjioned  im¬ 
mediately  to  the  ablating  surface  of  the  composite.  Therefore,  at  t  >5  s  a  profile 
of  tangential  stresses  has  two  domains  of  negative  values  (where  tempera¬ 

tures  are  relatively  low;  295K  <  6  <  1500A"  and  are  very  high  6  >  2400  K)  and 
one  domain  of  positive  values  (where  temperature  is  equal  to  1500  K  <  0  <2400 
K).  The  appearance  of  thermal  stresses  of  such  type  is  a  specific  featiure  of  the 
considered  class  of  Carbon/Phenolic  Composite. 

The  most  dangerous  stresses,  from  the  point  of  view  of  destruction,  are  tensile 
stresses  The  degree  of  danger  is  characterized  by  the  damage  parameter  Z2‘ 
the  closer  its  value  is  to  1,  the  closer  stress  is  to  ultimate  strength  a  of  the 
composite  FRFC  in  tension  under  the  given  tenperature  0.  Distributions  of  the 
damage  parameter  2:2 (r,  t)  versus  the  shell  thickness  for  different  times  t  axe  shown 
in  Figure  26.  Functions  Z2{r^  t)  have  a  typical  peak  at  t  >2  s,  which  moves  towards 
the  internal  surface  of  the  cylinder,  while  the  composite  is  warmed  up.  Value  of 


Distribution  of  volumetric  content  of  polymer  phase  ip2 
versus  thickness  of  the  cylindrical  shell  of  BRFC  at  different  times  t 
(r  -  radius  of  the  shell).  Symbol  *—  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Rnit)  at  times  t. 


Z2  reaches  its  maximum  at  the  top  of  this  peak  and  becomes  equal  to  0.40  -  0.75. 
And  absolute  maximum  of  value  Z2  equal  0.75  is  at  time  t  «  1  s,  when  minimal 
values  of  £22^  are  realized  at  the  ablating  surface  (shrinkage  reaches  its  maximum 
value).  Apparently,  this  time  moment  for  the  considered  type  of  BRFC  is  the  most 
dangerous  from  the  point  of  view  of  possible  destruction  of  the  shell  by  the  type 
(A). 

Besides  tangential  stresses,  radial  thermostresses  Cr  prove  to  be  also  very 
essential  in  the  cylindrical  shell  of  BRFC.  Fihure  27  shows  distributions  of  stresses 
ar{r,  t)  versus  the  shell  thickness  for  different  times  t.  During  initial  time  period 
of  heating  {t  <  0.2  s),  when  the  process  of  volumetric  ablation  has  no  time  to  be 
realized,  heat  deformation  ln^  is  positive  and  stresses  <7r  are  also  positive  (i.e. 
tensile).  The  process  of  volumetric  ablation  near  the  heated  surface  at  t  >  0.5 
s  leads  to  the  appearance  of  shrinkage  deformations  <  0)  and  compressive 

stresses  (Tr  within  the  shell  (Fig.  27).  However,  pore  gas  pressure  p  also  has  an 
essential  effect  on  a  character  of  radial  stresses  (Tr.  Within  the  zone  of  the  heated 
composite  (at  0  >  600  K),  as  it  was  mentioned  above,  there  appears  a  peak  of 
pore  pressure.  This  peak  causes  the  appearance  of  a  peak  (local  maximum)  of 
radial  stresses  <7^.  Value  of  stresses  (Tr  tbe  top  of  the  peak  can  be  both  positive 
and  negative.  As  seen  from  Fig.27,  at  times  t  =  0.5  and  7  s  maximum  of  stresses 
(Tr  is  positive,  i.e.  are  tensile  within  the  domain  ((Tr  ^  0.4  MPa),  and  at  time 
1  <  t  <  7  s  the  (Tr  stress  peak  is  negative  (min  (Tr  ^  -1.5  MPa),  and  the  (Tr  stress 
peak  is  displaced  towards  the  internal  surface,  while  the  composite  is  warmed  up. 

Danger  of  stresses  from  the  point  of  view  of  possible  destruction  of  com¬ 
posite  BRFC  by  the  type  (C)  (delamination)  is  described  by  damage  parameter  zi. 
Distributions  of  damage  parameter  zi  (r,  t)  versus  the  shell  thickness  for  different 
times  are  shown  in  Figure  28.  As  seen  from  this  figure,  a  character  of  distribution 
of  parameter  zi  has  a  peak  form.  An  appearance  of  peaks  is  caused  by  the  presence 
of  peaks  of  pore  gas  pressure  p.  As  follows  from  the  calculated  results,  maximum 
value  of  zi  is  very  high:  0.5  -  0.95;  it  means  that  a  state  of  composite  BRFC  is 
close  to  destruction  by  type  (C)  (due  to  delamination). 

Let  us  note  another  peculiarity.  For  composites  BRFC  with  ^  >  0°  for  non¬ 
zero  stresses  ar  there  appear  not  only  normal  stresses  an  but  also  stresses  of 
interlayer  shear  a^s,  which  also  raise  a  danger  of  delamination  of  the  composite. 


Distribution  of  tangential  stresses  oq  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t 
(r  -  radius  of  the  shell).  Symbol  <—  shows  a  location 
of  the  carried  off  surface  of  the  shell  ra  =  Roit)  at  times  t. 
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Fig.  26. 

Distribution  of  damage  parameter  Z2  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  / 

(r  -  radius  of  the  shell).  Symbol  <—  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Roit)  times  t. 


Fig.  27. 

Distribution  of  stresses  cr^  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t 
(r  -  radius  of  the  shell). 
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Fig.  28. 

Distribution  of  damage  parameter  z\  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t 
(r  -  radius  of  the  shell).  Symbol  shows  a  location 
of  the  carried  off  surface  of  the  shell  =  Roi^)  at  times  ti 
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It  is  due  to  this  that  in  Figure  28  function  Z2{r,t)  has  a  ’plateau’  on  the  right  side 
of  the  main  peak  of  values,  this  plateau  is  of  a  sufEciently  high  level  of  values; 
zi  =  0.20  -  0.35. 

Values  of  parameter  zi  reaches  their  maximum  near  the  ablating  surface  r  = 
Rjj  at  t  =0.5  s  and  at  final  time  of  calculation  t  =  7  s  near  the  internal  surface 
of  the  shell  r  =  .  Hence,  the  appearance  of  delaminations  within  a  structure  of 

the  considered  type  is  possible  near  both  the  external  and  internal  surfaces  of  the 
shell. 


10.4.2.  THE  EFFECT  OF  ANGLE  V’ 

Let  us  investigate  now  the  effect  of  angle  V’  on  a  thermomechanical  behaviour 
of  the  shell  of  BRFC  under  the  action  of  a  high-temperature  gas  flow. 

Calculations  being  similar  to  the  above-described  ones  were  conducted  for 
a  composite  with  angle  V’  =  0°,  i.e.  for  2-D  fabric  composite,  layers  of  which 
are  laid  perpendicularly  to  the  normal  vector  to  the  lateral  surface  of  the  shell. 
Calculations  in  Figs.  29  -  32  correspond  to  this  case. 

Figure  29  shows  distributions  of  temperature  9{r,  t)  versus  the  shell  thickness 
under  the  same  conditions  of  heating  as  in  Fig.22.  As  seen  from  the  figure,  warming 
up  the  shell  for  case  V’  =  0°  is  lower:  at  f  =  7  s  temperature  of  the  internal  surface 
does  not  practically  change  dg  =  295  K.  It  is  connected  with  that  heat  conductivity 
of  composite  ku  for  V’  =  0°  is  less  than  kn^  at  ^  =  30°  (see  Figures  2,  3). 

Lower  warming  up  the  composite  maJces  it  possible  that  gases  generating  in 
pores  at  volumetric  ablation  have  time  to  filtrate  and  outcome  partially  into  the 
surroundings.  Therefore  the  level  of  pore  pressure  (see  Fig.  30)  for  ^  =  0°  is  lower 
than  for  ^  =  30°.  Maximum  values  of  p  are  equal  to  45  and  100  atm,  respectively. 

Lowering  the  level  of  pore  pressure  p  leads  to  decreasing  peaks  of  radial  stresses 
(Fig.  31)  and  damage  parameter  zi  (Fig.  32).  Maximum  values  of  parameter  z\ 
for  V*  =  0°  does  not  exceed  0.1,  i.e.  a  danger  of  the  appearance  of  delaminations 
in  a  composite  with  ^  =  0°  proves  to  be  lower,  in  spite  of  higher  strength  of 
composite  BRFC  with  tp  =  30°. 

At  the  same  time  the  effect  of  angle  ^  on  a  character  of  tangential  stresses 
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<7^  and  damage  parameter  Z2  is  not  essential  practically;  maximum  level  of  values 
of  parameter  Z2  is  also  equal  0.75.  Hence,  the  main  cause  of  possible  destruction 
of  2-D  composite  with  V’  =  0°  is  cracking  of  fibres  near  the  ablating  surface  (type 

(A)). 

It  should  be  mentioned,  that  the  composite  with  ■0  =  30®  has  a  less  linear 
rate  of  recession:  D^i  =  0.57  mm/s  as  compared  with  D  =  0.61  mm/s  for  0  =  0° 
(see  Table  23)  and  also  a  less  mass  rate  of  recession  m  (see  Figs.  33  and  34)  and 
less  dimensionless  rate  of  recession  G.  Mass  rate  m  is  determined  as  follows: 

dAm 


where 


Am{t)  =  mo  —  m(<). 


is  a  loss  of  specific  mass  of  the  composite,  where  mo  and  m{t)  axe  determined  as 
follows: 

nT-o  =  Pq{R2  —  7?i); 

•/?£>(<) 


m{t)  = 


1  n 


p(r)rdr. 


here  p  —  Pi^fi  is  the  density  of  the  composite. 

Dimensionless  rate  G  is  determined  as  follows: 


G  = 


m 


(alcpY 


In  spite  of  that  ablation  rate  (linear  ,  mass  m  and  dimensionless  G)  of  com¬ 
posite  BRFC  with  0  =  30°  is  less  than  for  angle  0  =  0°,  temperatiure  of  the 
internal  surface  9s  proves  to  be  higher  (at  <  =  7  s  it  is  equal  to  414  and  293  K, 
respectively).  There  appears  a  natural  question:  what  variant  is  more  effective. 

If  we  estimate  these  variants  with  the  help  of  well-known  criterion  -  effective 
enthalpy  Jg//  determined  as 


then  we  obtain 


leff  —  12.27  MJ/kg  at  0  =  0°; 
Jg// =  12.95  MJ/kg  at  0  =  30°; 
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Fig.  30. 

Distribution  of  pore  gas  pressure  p  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t  and  ^  =  0° 
(r  -  radius  of  the  shell).  Symbol  <—  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Roit)  S'f  times  t'. 
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i.e.  variant  with  V’  =  30°  would  be  considered  more  eflfective.  However,  this 
criterion  does  not  take  into  account  that  temperature  6^  is  higher  for  rj)  =  30°.  In 
fact,  it  is  a  deficiency  of  the  criterion  of  effective  enthalpy. 

More  convenient  criterion  for  such  situations  is  another  one:  the  criterion  of 
dimensionless  parameter  of  efficiency  Di  introduced  in  [5]: 


Di  = 


Pqc{R2  -  RiYSa 
In  accordance  with  this  criterion: 


Di  =  15.5  for  V»  =  0“; 

Di  =  11.0  for  V’  =  30°; 

i.e.  composite  BRFC  with  ^  =  0°  proves  to  be  more  effective. 

If  one  have  no  requests  to  composite  BRFC  differing  from  requests  consid¬ 
ered  in  the  present  work,  then  the  criterion  on  efficiency  Di  is  more  preferable  as 
compared  with  leffi  and,  hence,  BRFC  with  V*  =  0°  should  be  considered  as  more 
effective. 
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Fig.  31. 

Distribution  of  stress  <Tr  versus  thickness 
of  the  cylindrical  shell  of  BRFC  at  different  times  t  and  ^  =  0° 
(r  -  radius  of  the  shell).  Symbol  <—  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Rnit)  times  t. 


Distribution  of  damage  parameter  zi  versus  thickness 
the  cylindrical  shell  of  BRFC  at  different  times  t  and  tp  =  0 
(r  -  radius  of  the  shell).  Symbol  shows  a  location 
of  the  carried  off  surface  of  the  shell  r2  =  Roit)  at  times  t. 
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Fig.  33. 

Dependence  of  mass  loss  Am  of  BRFC 
on  time  t  of  heating. 
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Varying  of  Loss  mass.  Mass  Ablatiom  Rate  and 
Dimensionless  Mass  Ablation  Rate  on  Heating  Time  t 
of  BRFC  with  ^=0° 


Time  t,s 

Loss  of 

Mass  Ablation 

Dimensionless 

Rate  m 

Mass 

Ablation 

kg/m2 

kg/(m2's) 

Rate 

G 

.1 

.0 

.000 

.000 

.2 

.1 

.  783 

.356 

.3 

.2 

.977 

.444 

.4 

.3 

.990 

.450 

.5 

.4 

1.061 

.483 

.6 

.5 

1.093 

.497 

.  7 

.6 

1.093 

.497 

.8 

.8 

1.093 

.497 

.9 

.9 

1.069 

.486 

1.0 

1.0 

1.089 

.495 

1.2 

1.2 

1.080 

.491 

1.4 

1.4 

1.067 

.485 

1.6 

1.6 

1.050 

.477 

1.8 

1.8 

1.048 

.476 

2.0 

2.0 

1.048 

.476 

2.2 

2.2 

1.037 

.471 

2.4 

2.4 

1.026 

.469 

2.6 

2.7 

1.023 

.466 

2.8 

2.9 

1.023 

.465 

3.0 

3.1 

1.016 

.462 

3.2 

3.3 

.998 

.458 

3.4 

3.5 

.997 

.454 

3.6 

3.7 

.997 

.453 

3.8 

3.9 

.991 

.450 

4.0 

4.1 

.987 

.449 

4.2 

4.3 

.983 

.447 

4.4 

4.4 

.981 

.446 

4.6 

4.6 

.970 

.441 

4.8 

4.8 

.969 

.441 

5.0 

5.  0 

.963 

.438 

5.2 

5.2 

.954 

.434 

5.4 

5.4 

.954 

.434 

5.  6 

5.6 

.946 

.430 

5.8 

5.8 

.942 

.429 

6.0 

6.0 

.942 

.428 

6.2 

6.2 

.928 

.422 

6.4 

6.4 

.928 

.422 

6.6 

6.5 

.928 

.422 

6.8 

6.  7 

.927 

.421 
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Varying  of  Loss  mass.  Mass  Ablatiom  Rats  and 
Dimensionless  Mass  Ablation  Rate  on  Heating  Time  t 
of  BRFC  with  T=30° 


Time  t,s  Loss  of  Mass  Ablation  Dimensionless 

mass  -Am,  Rate  m  “  Mass  Ablation 

kg/m2  kg/(m2-s)  Rate  G 


.1 

.0 

.2 

.1 

.3 

.3 

.4 

.4 

.5 

.5 

.6 

.6 

.  7 

.  7 

.8 

.9 

.9 

1.0 

1.0 

1.1 

1.2 

1.3 

1.4 

1.6 

1.6 

1.8 

1.8 

2.0 

2.  0 

2.3 

2.2 

2.5 

2.4 

2.  7 

2.6 

2.9 

2.8 

3.2 

3.0 

3.4 

3.2 

3.6 

3.4 

3.8 

3.6 

4.  0 

3.8 

4.2 

4.0 

4.4 

4.2 

4.6 

4.4 

4.9 

4.6 

5.1 

4.8 

5.3 

5.0 

5.5 

5.2 

5.  7 

5.4 

5.9 

5.6 

6.1 

5.8 

6.3 

6.0 

6.5 

6.2 

6.7 

6.4 

6.9 

6.6 

7.1 

6.8 

7.3 

7.0 

7.4 

.  000 

.000 

.  886 

.403 

1.207 

.549 

1.221 

.555 

1.231 

.560 

1.220 

.555 

1.216 

.553 

1.216 

.553 

1.196 

.544 

1.201 

.546 

1.183 

.538 

1.169 

.531 

1.155 

.525 

1.143 

.520 

1.129 

.513 

1.113 

.506 

1.107 

.503 

1.103 

.502 

1.099 

.500 

1.088 

.494 

1.075 

.488 

1.070 

.487 

1.068 

.486 

1.059 

.481 

1.046 

.475 

1.046 

.475 

1.042 

.474 

1.031 

.468 

1.025 

.466 

1.023 

.465 

1.015 

.461 

1.008 

.458 

1.006 

.457 

.997 

.453 

.996 

.453 

.990 

.450 

.983 

.447 

.982 

.446 

.977 

.444 

.978 

.444 
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Varying  of  Recession  Rate,  Temperature  of  Ablation  Surface 
and  Temperature  of  Internal  Surface  vs  Heating  Time  t 
of  BRFC  with  F'O® 


Time  t,s 

Rate  of 

Temperature 

Temperature 

Recession 

of  Ablation 

of  Internal 

mm/s 

Surface 

Surface 

.1 

.000 

2142 

294 

.2 

.363 

2911 

294 

.3 

.558 

2972 

294 

.4 

.574 

2988 

294 

.5 

.579 

3004 

294 

.6 

.583 

3014 

294 

.  7 

.586 

3024 

294 

.8 

.589 

3030 

294 

,9 

.591 

3036 

294 

1.0 

.593 

3041 

294 

1.2 

.595 

3049 

294 

1.4 

.598 

3055 

294 

1.6 

.599 

3060 

294 

1.8 

.601 

3066 

294 

2.0 

.602 

3070 

294 

2.2 

.604 

3073 

294 

2.4 

.604 

3075 

294 

2.6 

.605 

3078 

294 

2.8 

.606 

3080 

294 

3.0 

.607 

3082 

294 

3.2 

.607 

3084 

294 

3.4 

.608 

3085 

294 

3.6 

.608 

3087 

294 

3.8 

.609 

3088 

294 

4.0 

.609 

3092 

294 

4.2 

.610 

3092 

294 

4.4 

.610 

3093 

294 

4.  6 

.611 

3094 

294 

4.8 

.611 

3095 

294 

5.0 

.611 

3096 

294 

5.2 

.612 

3097 

294 

5.4 

.612 

3098 

294 

5.  6 

.612 

3099 

294 

5.8 

.612 

3099 

294 

6.  0 

.613 

3100 

294 

6.2 

.613 

3101 

294 

6.4 

.613 

3101 

294 

6.  6 

.613 

3102 

294 

6.8 

.613 

3102 

295 

7.0 

.613 

3102 

295 
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Varying  of  Recession  Rate,  Temperature  of  Ablation  Surface 
and  Temperature  of  Internal  Surface  vs  Heating  Time  t 
of  BRFC  withy^SO'^ 


Time  t,s 

Rate  of 

Temperature 

Temperature 

Recession 

of  Ablation 

of  Internal 

Difr^  mm/s 

Surface  i^,K 

Surface  ,K 

.1 

.000 

2142 

293 

.2 

.344 

2911 

293 

.3 

.520 

2972 

293 

.4 

.536 

2988 

293 

.5 

.540 

3004 

293 

.6 

.544 

3014 

293 

.7 

.547 

3024 

293 

.8 

.549 

3030 

293 

.9 

.551 

3036 

293 

1.0 

.553 

3041 

293 

1.2 

.555 

3049 

293 

1.4 

.557 

3055 

293 

1.6 

.559 

3060 

293 

1.8 

.561 

3066 

293 

2.0 

.562 

3070 

293 

2.2 

.563 

3073 

293 

2.4 

.563 

3075 

293 

2.6 

.564 

3078 

293 

2.8 

.565 

3080 

293 

3.0 

.565 

3082 

293 

3.2 

.566 

3084 

293 

3.4 

.566 

3085 

293 

3.6 

.567 

3087 

293 

3.8 

.567 

3088 

293 

4.0 

.568 

3092 

293 

4.2 

.569 

3092 

293 

4.4 

.569 

3093 

293 

4.6 

.569 

3094 

294 

4.8 

.570 

3095 

294 

5.0 

.570 

3096 

296 

5.2 

.570 

3097 

297 

5.4 

.570 

3098 

300 

5.6 

.571 

3099 

304 

5.8 

.571 

3099 

310 

6.0 

.571 

3100 

318 

6.2 

.571 

3101 

328 

6.4 

.571 

3101 

342 

6.  6 

.572 

3102 

361 

6.  8 

.572 

3102 

384 

7.0 

.572 

3103 

414 

Chapter  XI 


THE  EFFECT  OF  INCLINATION  ANGLE  ON 
PROPERTIES  OF  CARBON/PHENOLIC  BRFC 


Table  25  shows  summary  results  of  all  calculations  conducted  for  the  present 
report  for  properties  of  Carbon/Phenolic  BRFC  at  two  angles  =  0  and  30°.  The 
purpose  of  the  table  is  to  investigate  the  effect  of  inclination  angle  ^  on  properties 
of  composites  and  on  a  thermo-mechanical  behaviour  of  a  cylindrical  shell  made 
of  BRFC  under  external  heating. 

In  the  last  column  of  the  table  there  are  relative  changes  (m  %)  of  values  of 
the  functions  at  passage  from  angle  =  0°  to  t/>  =  30°.  Sign  'pins'  means  an 
increase  of  values,  ’minus’  means  a  decrease. 
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temperature 

yr  —  Qo  y- 

1 .  Across  Ply  Thermal 
Conductivi  ty  {  W/m  K) 

295  K 
3000  K 

.65 

5.  72 

1.10 

10.40 

+69 

+82 

2.  Longitudinal  Thermal 
Conductivity  (Vl/m  K) 

295  K 
3000  K 

2.25 

25.00 

2.00 

19.  76 

-11 

-21 

3.  Across  Ply  Thermal 
Expansion  /  % 

600  K 
3000  K 

1.30 

-10.60 

.976 

-8.05 

-25 

-20 

4 .  Longi tudinal  Thermal 
Expansion  ,  % 

295  K 
3000  K 

.01 

.14 

.014 

-2.568 

+40 

5.  Across  Ply  Elastic 
Modulus  ,  GPa 

295  K 
3000  K 

6.  79 

0.23 

S.34 

0.f4 

-Xi 

-t4 

6.  Longitudinal  Elastic 
Modulus  ,  GPa 

295  K 
3000  K 

60.23 

30.00 

4.41 

O.  22 

-99 

7.  Shear  Elastic 

Modulus  ,  GPa 

295  K 
3000  K 

1.679 

.083 

4.189 

.200 

+150 

+141 

8.  Strength  in  Tension 
along OXi  ,  ^5*^^  ,  MPa 

295  K 
3000  K 

11.58 

.46 

14.92 

.60 

+29 

+30 

9.  Strength  in  Tension 
along  Oxj  , 

295  K 
3000  K 

399. 70 

76.  80 

23.20 

.90 

-94 

-99 

10.  Shear  Strength 

,  MPa 

295  K 
3000  K 

25.14 

1.60 

12.92 

.52 

-49 

-55 

11.  Linear  Ablation 

Rate 

3000  K 

Pa =  3  atm 

(  ^/Cf,  )  =2. 

.61 

2 

.57 

-6.5 

12.  Temperature  ^s-^o  (K) 
of  internal  'cold'  surface 
after  fc=:  7  s  of  heating 

12  MW/m  0 

k-  12  mm 

119 

13.  Temperature  &uf  of 
ablation  surface,  K 

=  12  MW/m  3102 

Pa  -  3  atm 

(  ^  /Cf,)  =2.2 

3102 

14.  Mass  ablation  rate 
/n  ,  kg/(m  ^  s) 

-  ff  - 

.978 

.927 

-5.5 

15.  Dimensionless_mass 
ablation  rate  & 

-  ff  - 

.44 

.42 

-4.8 

16.  Maximal  value  of  damage 
parameter  Z^in  cylindrical  , 

-  fr  . 

shell 

.1 

.95 

+850 

17.  Maximal  value  of  damage 
parameter  in  cylindrical  ; 

.  ff  . 

shell 

.  75 

.  75 

0 

18.  Maximal  value  of  radial 
thermal  stresses  in  cylindr 
ical  shell  ,  MPa 

.  n  . 

-0.25 

.4 

19.  Maximal  tangential  thermal  -  »  - 
stresses  in  cylindrical  shell 
,  MPa 

57 

57 

0 

20.  Maximal  value  of  pore 
pressure  p,  MPa 

4.5 

10 

+122 

21.  Effective  enthalpy Jg^,  MJ/kg  -  "  ■ 

12.27 

12.95 

■tS.S 

22.  Dimensionless  parameter 

-  n  - 

15.5 

11. 
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CONCLUSIONS 


1.  In  the  present  Report  the  model  developed  before  on  thermo- mechanical  processes  of 
ablation  for  Carbon/Phenolic  Bias- Reinforced  Composites  under  high  temperatures 
has  been  modified. 

2.  Calculations  have  been  conducted  for  heat-physical  characteristics  of  BRFC  (heat 
conductivity,  heat  expansion  (shrinkage))  for  different  inclination  angles  as  initial 
data  we  have  used  experimental  results  obtained  in  [11]. 

It  is  established  that  for  V*  =  30°  Across  Ply  Conductivity  is  higher  by  70  -  80  % 
than  for  BRFC  with  ^  =  0°.  Due  to  this,  temperature  on  the  internal  wall  of  the 
cylindrical  shell  with  angle  ip  =  30°  under  external  heating  by  flux  =  12  MW/m^ 
will  be  higher  {$3  —  OofnO  and  119  °C,  respectively,  at  time  t  =  7  s  of  heating). 

It  is  also  established  that  Longitudinal  Thermal  Expansion  l22t/>  being  the  cause 
of  appearance  of  tangential  stresses  in  the  cylindrical  shell  for  ^  =  0  and  30°  is 
practically  the  same. 

At  the  same  time  heat  deformation  £33^  essentially  increases  in  its  absolute  magni¬ 
tude,  but  its  role  for  the  problem  on  heating  of  a  cylindrical  shell  is  not  considerable. 

3.  With  the  help  of  the  modified  model,  calculations  have  been  conducted  for  thermo¬ 
mechanical  characteristics  (elastic  moduli,  strength)  of  BRFC  for  different  incli¬ 
nation  angles  ip.  It  has  been  established  that  for  values  ip  =0  and  30°  Across  Ply 
Elastic  Modulus  and  Longitudinal  Elastic  Modulus  E2,/,  are  practically  the  same 
within  all  the  temperature  interval,  only  longitudinal  modulus  E3  decreases  that  has 
no  essential  effect  on  the  problem  of  heating  a  cylindrical  shell. 

Increasing  angle  ip  from  0°  to  30°  leads  to  growing  Strength  in  Tension  along 
0x3  axis  by  30  %  and  Shear  Elastic  Modulus  by  150  %. 

4.  With  the  help  of  the  model,  numerical  calculations  of  linear  rate  of  ablation  of 
BRFC  have  been  conducted  for  different  angles  ip.  Comparison  with  experimental 
values  taken  from  scientific  literature  data  for  ip  =  0°  has  shown  a  good  accuracy  of 
prediction  of  the  linear  ablation  rate  with  the  help  of  the  model.  It  has  been  set  that 
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a  contribution  of  thermo-mechanical  erosion  into  the  total  linear  rate  of  ablation  is 
very  essential:  for  considered  conditions  of  heating  (9®  =  12  MW/m^,  Po  =  3  atm) 
it  is  equal  to  50  -  70  %. 

It  has  been  also  established  that  an  increase  of  angle  tp  from  0°  to  30°  leads  to 
decreasing  the  summarized  linear  rate  of  ablation  by  7  %,  the  mass  rate  of 
ablation  m  by  6  %,  dimensionless  rate  of  recession  G  by  5  %;  while  temperature  of 
the  heated  surface  9yj  does  not  change  practically. 

5.  With  the  help  of  the  model,  calculations  have  been  conducted  for  thermo-stresses 
and  internal  heat  and  mass  transfer  in  a  cylindrical  shellmade  of  BRFC  with  angles 
V’  =  0°  and  30°  under  external  convective  heating  with  parameters  =  12  MW /m^ 
and  Po  =  3  atm. 

It  has  been  established  that  temperature  on  the  internal  nonheated  surface  of  the 
composite  grows  with  increasing  angle  ip,  and  maximum  of  pore  gas  pressure  p  also 
increases  (more  by  two  times). 

6.  Increasing  angle  ip  has  no  effect  practically  on  a  level  of  tensile  tangential  stresses 
(Tg  near  the  external  surface  of  the  cylindrical  shell:  their  maximum  is  equal  ap¬ 
proximately  to  57  MPa,  that  ensures  value  of  damage  parameter  equal  to  Z2  «  0.75. 
At  the  same  time  the  increase  of  angle  ip  from  0°  to  30°  leads  to  sharp  increasing 
radial  stresses  within  the  shell:  their  maximal  values  become  positive  (tensile) 
and  equal  to  «  0.4  MPa,  and  value  of  damage  parameter  z\  grows  up  to  0.95,  i.e. 
approaches  practically  the  ultimate  magnitudewhen  there  occurs  a  destruction  of 
the  composite  by  the  delamination  type. 

The  main  cause  of  growth  of  damage  parameter  zi  at  ip  =30°  is  higher  pore  gas  pres¬ 
sure  than  for  ip  =  0°  and  the  appearance  of  stresses  of  interlayer  shear  at  inclined 
location  of  the  fabric  layers  in  the  composite. 

7.  Thus,  the  mathematical  model  developed  of  coupled  thermo-mechanical  processes 
of  ablation  of  composites  under  the  action  of  high-temperature  gas  flows  allows  us: 

•  to  perform  calculations  of  heat-physical  and  thermo-mechanical  characteristics  of 
fabric  2-D  composites,  including  BRFC,  taking  as  initial  data  only  properties  of  the 
matrix,  fibres  and  parameters  of  their  geometrical  interlocaltion; 

•  to  conduct  calculations  of  ablation  rates  (mass  and  linear)  of  composites  under 
the  action  of  high-speed  hot  gas  flows,  including  account  of  both  thermo-chemical 
processes  (combustion,  sublimation,  pyrolysis)  and  thermo-mechanical  processes  (in¬ 
ternal  and  external  mechanical  erosion  of  different  types).  Here  for  calculations  of 


ablation  rates  of  composite  we  have  no  need  to  know  ablation  rates  of  its  matrix 
and  fibres:  they  can  be  determined  with  the  help  of  the  same  heat-physical  and 
mechanical  properties  of  the  matrix  and  fibres; 

•  to  conduct  calculations  of  thermo-stresses  and  parameters  of  internal  heat  and 
mass  transfer  in  structures  of  ablative  composites  under  the  action  of  gas  flows; 

•  to  optimize  an  internal  structure  and  composition  of  ablative  composites  in  order 
to  obtain  the  best  characteristics  and  the  highest  efficiency  of  the  materials. 

8.  Verification  of  the  developed  model  has  shown  that: 

•  the  model  ensures  a  quite  satisfactory  accuracy  of  prediction  of  the  characteristics 
of  ablative  composites:  the  best  estimations  are  at  room,  high  and  superhigh  tem¬ 
peratures,  the  worst  -  at  elevated  temperatures  (294  -  600  K);  the  best  prediction  is 
realized  for  elastic  moduli,  strength,  gas  permeability  and  heat  capacity,  the  worst 
-  for  heat  expansion  and  heat  conductivity,  it  is  explained  by  the  absence  of  suffi¬ 
ciently  reliable  data  on  thermal  expansion  and  conductivity  of  matrix  and  fibres; 

•  the  model  ensures  a  good  accuracy  of  prediction  of  linear  and  mass  rates  of  com¬ 
posite  ablation  under  the  action  of  gas  flows,  that  is  achieved  at  the  expense  of 
sufficiently  detailed  modelling  of  thermo-mechanical  erosion  processes  of  compos¬ 
ites  giving  a  main  contribution  into  the  ablation  process  at  great  values  of  velocity 
head  of  the  flow  p^. 

The  limited  scope  of  experimental  data  does  not  allow  us  to  conduct  more  detailed 
verification  of  the  model,  in  particular  for  different  values  of  velocity  head  Pq  and 
heat  flux  q°. 

9.  To  continue  the  work  in  future,  we  suggest  on  the  basis  of  the  present  model  to 
develop  a  method  of  optimal  designing  of  composition  and  internal  structure  of  ab¬ 
lative  composite  materials.  This  method  will  allow  us  to  select  such  composition 
of  composite  materials  (matrix,  one,  two,  three  or  more  fillers)  and  such  structure 
(angles  of  reinforcing,  volumetric  content  of  different  fillers  etc.)  which  would  en¬ 
sure  a  creation  of  ablative  composites  with  the  best  characteristics  and  the  highest 
efficiency  in  composition  of  structures  under  the  concrete  conditions  of  heating.  It 
should  be  noted  that  development  of  criteria  of  what  material  can  be  considered 
as  more  effective  will  be  the  most  important  problem  in  this  work.  For  example, 
under  the  certain  conditions  of  heating,  material  with  higher  heat  conductivity  and 
smaller  rate  of  ablation  can  prove  to  be  more  effective,  and  under  other  conditions 
of  heating  the  most  effective  material  can  have  opposite  properties. 


REFERENCES 


1.  ’Verification  of  Model  for  Heat-Mechanical  Processes  in  Ablating  Composites’. 
First  progress  report  F61708-96-W0227,  1996. 

2.  ’Development  of  Models  for  Heat-Mechanical  Processes  in  Thermal  Protection 
Made  of  Ablating  Composites’,  Report  No.l  on  the  contract  SPC-94-4104, 
1995. 

3.  ’Development  of  Models  for  Heat- Mechanical  Processes  in  Thermal  Protection 
Made  of  Ablating  Composites’,  Report  No.2  on  the  contract  SPC-94-4104, 
1995. 

4.  ’Development  of  Models  for  Heat-Mechanical  Processes  in  Thermal  Protection 
Made  of  Ablating  Composites’,  Report  No.3  on  the  contract  SPC-94-4104, 
1995. 

5.  Dimitrienko,  lon.I.,  ,  Efremov,  G.A.  and  Chernyavsky,  S.A.,  ’Optimum  De¬ 
sign  of  Erosion-Stable  Heatshield  Composite  Materials’.  Applied  Composite 
Materials,  v.4,  N  1,  pp. 35-52  (1997). 

6.  Dimitrienko  lou.I.  Thermal  stresses  and  heat-mass-transfer  in  ablating  com¬ 
posite  materials.  J.Heat  and  Mass  Transfer,  vol.38,  N  1,  1995,  pp. 139-146. 

7.  Dimitrienko  lou.L,  Effect  of  Finite  Deformations  on  Internal  Heat-Mass- 
Transfer  in  Elastomer  Ablative  Materials,  Int.  Journal  of  Heat  Mass  Transfer, 
vol.  40,  No.3,  1997,  pp. 699-709. 

8.  Dimitrienko  lou.I.,  Internal  Heat-Mass- Transfer  and  Stresses  in  Thin- Walled 
Structures  of  Ablating  Materials,  Int.  Journal  of  Heat  Mass  Transfer,  vol.40, 
No.7,  1997,  pp.1701-1711. 

9.  Dimitrienko  Yu.L,  Thermal  Stresses  in  Ablative  Composite  Thin- Walled 
Structures  under  Intensive  Heat  Flows,  Int.  Journal  of  Engineering  Science, 
vol.35,  N  1,  1997,  pp.15-31. 

10.  Dimitrienko,  lou.I., ,  Efremov,  G.A.,  Chernyavsky,  S.A.  and  Medvedev,  Yu.V., 
’Theory  and  Synthesis  of  Advanced  Thermal-Protective  Composite  Materials’, 
Applied  Composite  Materials,  v.  2,  N  6,  pp. 367-384  (1995). 

11.  ’Mechanical  and  Thermal  Data  Package  for  Generic  2-D  Carbon  Phenolic 
Composites.  Wright  Pat.  AFB.  Information  materials,  1996. 


12.  Dimitrienko  lou.L,  Efremov  G.A.,  Epifanovsky  I.S.  Reusable  Re-entry  Vehicles 
with  Reclaimable  Ablating  Thermal  Protection.  Preprints  of  19th  Int.Symp. 
on  Space  Technology  and  Science  May  1994,  Yokohama,  Japan,  94-b-26. 

13.  Polezhaev  Yu.V.,  Yurevich.  Thermal  Protection.  Moscow,  Energia,  1976  (in 
Russian). 

14.  Dimitrienko  Yu. I.,  Thermomechanical  Behaviour  of  Composite  Materials  and 
Structures  under  High  Temperatures.  1.  Materials.  Composites.  Part  A: 
Applied  Science  and  Manufacturing.  1997  (in  print). 

15.  Dimitrienko  Yu. I.,  Thermomechanical  Behaviour  of  Composite  Materials  and 
Structures  under  High  Temperatures.  2.  Structures.  Composites.  Part  A: 
Applied  Science  and  Manufacturing.  1997  (in  print). 

16.  Dimitrienko  Yu.I.,  Modelling  of  Mechanical  Properties  of  Composite  Materials 
under  High  Temperatures.  Part  1.  Matrix  and  Fibres.  Applied  Composite 
Materials.,  1997  (in  print). 

17.  Dimitrienko  Yu.I.,  Modelling  of  Mechanical  Properties  of  Composite  Materi¬ 
als  under  High  Temperatures.  Part  2.  Unidirectional  Composites.  Applied 
Composite  Materials.,  1997  (in  print). 

18.  ’Engineering  Design  for  Plastics’.  Ed.  by  E.Bacr,  New  York,  Reinhold  Publ. 
Corp.,  1967. 

19.  Dimitrienko  Yu.I.,  Modelling  of  Erosion  Combustion  of  Energetic  Materials  in 
High-Enthalpy  Flows,  Combustion  And  Flame,  1997  (in  print). 

20.  Laub  B.,  Thermo-Chemical  Ablation  of  Tantalum  Carbide  Loaded  Carbon- 
Carbons.  AIAA  Pap.,  1980,  N  1476. 

21.  Pobedrya  B.E.,  Lectures  on  Tensor  Analysis.  Moscow  State  University  Publ., 
1987  (in  Russian). 


